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We obtain analytic expressions for the time correlation functions of a liquid of spherical particles, 
exact in the limit of high dimensions d. The derivation is long but straightforward: a dynamic virial 
expansion for which only the first two terms survive, followed by a change to generalized spherical 
coordinates in the dynamic variables leading to saddle-point evaluation of integrals for large d. The 
problem is thus mapped onto a one-dimensional diffusion in a perturbed harmonic potential with 
colored noise. At high density, an ergodicity-breaking glass transition is found. In this regime, our 
results agree with thermodynamics, consistently with the general Random First Order Transition 
scenario. The glass transition density is higher than the best known lower bound for hard sphere 
packings in large d. Because our calculation is, if not rigorous, elementary, an improvement in the 
bound for sphere packings in large dimensions is at hand. 


Introduction - The physics of liquids and glasses be¬ 
longs to the group of fields that are victims of the lack 
of a small parameter. Many approximations have been 
proposed over the years, but they suffer from the uncer¬ 
tainty about what is the limit in which they are sup¬ 
posed to become exact. This has been true both for 
equilibrium and for dynamic properties. From the point 
of view of dynamics, an extreme case is that of Mode- 
Coupling Theory (MCT) [1-3]: it may be introduced 
by an (uncontrolled) resummation of an infinite subset 
of diagrams. The Mode-Coupling approximation yields 
Mode-Coupling dynamics: the phenomenology depends 
on the approximation itself [4], somewhat like a harmonic 
approximation is expected to predict oscillations. 

An often used remedy to the absence of a small pa¬ 
rameter [5] is to promote the system to d dimensions, 
solve the large d limit, and (eventually) expand around. 
This strategy has been used with success for liquids [6-8], 
strongly coupled electrons [9], atomic physics [10], gauge 
field theory [11], and most recently the thermodynam¬ 
ics of amorphous systems [12, 13]. In this paper we 
extend this procedure to the dynamics of liquids made 
of spherical particles. We restrict ourselves to equilib¬ 
rium, although the extension to the glassy off-equilibrium 
“aging” regime is at hand. It has been a long-standing 
question whether MCT becomes exact in infinite dimen¬ 
sion [4, 14-18], and the present computation gives an 
answer. 

Statement of the main result - We shall consider a 
system of N identical particles, interacting via a spher¬ 
ical potential V(r) of typical interaction length cr in c? 
dimensions, and obtain a solution for the equilibrium 
time correlations of the resulting liquid, that becomes 
exact in the limit d —i (x. We need to confine the par¬ 
ticles in a finite volume V. It is very convenient to do 
this in such a way that the “box” does not break ro¬ 
tational and translational invariances, which are crucial 
in our developments. A practical way to do this is to 
consider particles living on points Xi (i = on 

the d-dimensional surface of a hypersphere with Xi- Xi = 
= cr2Aiiq/(2d) > cr^ (/X = 1, • • • , d -b I). 
The thermodynamic limit i? —)■ oo with constant density 


p = N/V, in which the flat space is recovered, will be 
taken before d —> oo. Rotations and translations in d di¬ 
mensions - with dimensions d(d—1)/2 and d, respectively 
- are encoded in the rotations in d -I- 1 dimensions, with 
dimension d{d+ l)/2. We consider a Langevin dynamics 


mxi{f) + 7ii(t) = -u^{f)xi{t) - Vjqid -b f^{f) , (I) 

where is a white noise with {ffit)fj {t')) = 

2Tj6ijSf_,^>6{t-t'), H = and T = l/jS is 

the temperature. Here, and in what follows, (•) denotes 
average over noise fi and/or initial conditions. The Vi are 
Lagrange multipliers, imposing the spherical constraints 
Xi ■ Xi = 0. For d —>■ oo, they do not fluctuate and their 
value Vi ^ V = dTp/R^ is proportional to the equilib¬ 
rium reduced pressure p = fdP/p [19]. We shall in what 
follows treat the overdamped case m = 0, but the inertial 
term may be reinstalled at any stage (and in that case 
the thermal bath may be disconnected setting 7 = 0, to 
recover a purely Newtonian dynamics). 

We define the adimensional scaled correlation and re¬ 
sponse functions (see Fig.l in the Appendix): 


C{t,t') ^ C,{t,t') = 
R{t,t') ~ R^{t,t') = 

( 7 ^ 


-ix.it) ■ x,it') , 

Cr2 ^ dhfit') ’ 

Xi{t')\'^ = Aliq - C{t,t') , 


( 2 ) 


where hi is an external field conjugated to Xi, and we 
have also introduced the single particle mean square dis¬ 
placement A [19]. As we shall see, these quantities have 
a finite and non-fluctuating limit for d —)■ 00. In the 
following and in the Appendix, we shall derive exact 
equations for the correlations (2). For simplicity, we will 
restrict to equilibrium where C{t,t') = C{t — t'), and 
R{t) = —l39{t)C{t); here 9{t) is the Heaviside step func¬ 
tion. 

Our main result is that correlation functions are ob¬ 
tained in terms of a memory kernel M through the fol- 
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lowing equations: 

^C{t) = ~ ~ u)C(u) , 

J, ^ ( 3 ) 

^A{t) = T — P du M{t — u)A{u) , 

Jo 

the second relation being valid for A ^ Auq. 

To give the expression of the memory kernel M, let us 
define a scaled density = pVdicr)/d where Vdiu) is the 
volume of a d-dimensional sphere of radius u, a scaled 
7 = a scaled potential V{y) = V{<j{l-\-y/d)), 

and the inter-particle force F{y) = —V'{y). M is then 
obtained in terms of a one-dimensional effective dynamics 
with a colored noise 

jy{t) = -w'{y{t)) - P du M{t - u)y{u) + C{t) , 

(C(t)C(f')) = - t') + Mit - t') , 

from the self-consistent equation: 

M{t - i') = f / dyo e-^-^y-HF{t)F{t'))M,yo ■ (5) 

The average in (5) is over the process (4) with y{t = 0) = 
7 /0- The total effective potential has the form (see Fig.4 
in the Appendix): 

Tt/^ 

w{yo) = V{yo)-Tyo +. (6) 

The quadratic part of the potential plays the role of the 
confining “box”. In fact, it is negligible for finite times 
and large Aiiq, and the probability distribution is ex¬ 
ponential in y near j/ ~ 0 (expressing the growth of en¬ 
tropy as a function of distance along the d-f 1-dimensional 
sphere), the region relevant for those times. Reassuringly, 
box details are irrelevant at short times. 

Finally, Eq. (5) has a microscopic counterpart (cf. 
Appendix); in fact, the function M{t) is, in the large 
d limit, the autocorrelation of the inter-particle forces 

Fijit) = -'^V{Xi{t) - Xj{t)), 

= , ( 7 ) 

i¥=3 

a result that provides a physical interpretation of M{t). 

Relation with MCT - If M were a simple function of 
C, M = F(C'), then Eq. (3) would be in the schematic 
MCT form. As is well known, schematic MCT is ob¬ 
tained as the exact dynamics of a system of spherical 
spins sf = N with p-spin random interactions [20- 
22 ], for which 

M = F(C) = . (8) 

However, as soon as one considers non-spherical vari¬ 
ables, e.g. soft-spins with a potential 14(s) = a(s^ — 1)^, 


one obtains an equation like (4) with this I4(s) [21, 23, 
24]: 


s{t) =—V'(s) —/3 f du M{t — u)s{u) + C{t) . (9) 

Jo 

Here again, Eq. ( 8 ) holds and the system is closed by 
C{t — t') = {s{t)s{t')). Within the liquid phase, this 
more general form of dynamic equation has essentially 
the same phenomenology as schematic MCT. Our sys¬ 
tem of equations belongs to this more general class, and 
they thus show exactly the same MCT phenomenology 
for what concerns universal quantities that are indepen¬ 
dent of details of the memory kernel (e.g. the dynami¬ 
cal scaling forms and the relations between critical expo¬ 
nents) [3, 25, 26]. 

However, important quantitative differences are ob¬ 
served with respect to applying the MCT approximation 
to the intermediate scattering function, which leads to 
the “standard” formulation of MCT for liquids [1, 2]. 
Standard MCT has the same qualitative structure as 
schematic MCT, but also provides quantitative results 
for the self and collective scattering functions in all di¬ 
mensions, in particular in d = 3 [27, 28]; its d —)■ oo limit 
was discussed in Refs. [15, 16]. Our result in d —)■ oo 
is formulated in terms of A{t,t'), and most of the other 
natural observables are functionals of A{t, t'). For exam¬ 
ple, for qa ^ 1 , we have for the self intermediate 
scattering function (cf. Appendix): 


(l>\{t,t!) = exp 


9^2 
' 2d2 


A(t,C) 


( 10 ) 


in contrast to the non-Gaussianity in q one finds within 
MCT close to the plateau [15, 16]. One could then write 
our equations in terms of The result, however, is 
different from standard MCT and in particular our kernel 
M is not an analytic function of 

Because our equations fall in the same universality 
class as schematic MCT, but provide different quanti¬ 
tative results with respect to standard MCT in d —>■ oo, 
it remains a matter of taste if one wishes to call them 
with the same name or, more generally, “dynamic Ran¬ 
dom First Order Transition (RFOT)” [14]. 

Sketch of the derivation - Let us outline the main 
steps in the derivation; more details are given in the Ap¬ 
pendix. In order to construct the high-dimensional limit, 
a virial expansion is a reliable method. Following [29], we 
exploit the well-known analogy between trajectories and 
polymers. The dynamics are generated by a sum of tra¬ 
jectories in d-dimensional space, in our case, the surface 
of a d-|-1-dimensional sphere. To each trajectory is associ¬ 
ated an Onsager-Machlup probability weight which is the 
exponential of an action. The sum over all trajectories 
of this quantity is analogous to a partition function and 
is thus used to generate averages over the Langevin pro¬ 
cess (1). The action is expressed in terms of trajectories 
of the Xi and auxiliary “response” variables Xi (Martin- 
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Siggia-Rose-De Dominicis-Janssen generating path inte¬ 
gral) [20, 22]. The result reads, in the ltd convention: 


N 


i—1 

4.1x,il = 7/d((ri^ + ,-i.i + .x.x) , 

JT[a;i — X2, Xi — X 2 ] = i J dt (xi — X 2 ) ■ VR(|a;i — X2\) ■ 


( 11 ) 


Following standard liquid theory [19, 29], we introduce 
the density function for trajectories 

p[x,x] = (^^Y2d[x - x^]d[x - x^^ , (12) 


where (5[a;] is the functional Dirac S (a product of deltas 
over all times) and we construct a virial (Mayer) expan¬ 
sion as a power series in p[x,x]. One can show that all 
terms involving a product of more than two density fields 
are subleading for d —> 00 [8]: they are exponentially sup¬ 
pressed because of the requirement that the three tra¬ 
jectories overlap (see Fig.l in the Appendix), which is 
exponentially unlikely in large d. Truncating the virial 
expansion accordingly, we get (cf. Appendix): 


5 = 
-h 


In Zjv 
N 


N 

Y 


J Dp[xi,xi]Dp[x2,X2]f[xi - X2,Xi - X2] , 


(13) 


also a; = jxi — X 2 \'^, to = {xi — X 2 ) • {xi — X 2 ). The 
Jacobian may be calculated with the same methods (cf. 
Appendix), and the crucial result is that at the saddle- 
point level Qi = Q 2 = Q with Q determined by the 
same saddle-point as in Eq. (14), while the remaining 
integration over Lo{t),Ld(t) is effectively one-dimensional. 
Changing variables with uj{t) = ct^(1 J- y{t)/d) leads 
to a finite integration over y(t) that eventually gives 
Eq. (4). Hence, the typical distance between two trajec¬ 
tories turns out to be cr-1-0(1/d). This scaling physically 
tells us that a particle vibrate inside a cage with 1/d 
amplitude and interacts with 0(d) neighbors. Finally, 
equilibrium and causality at the saddle-point level im¬ 
ply D{t,t') = 0 and the Fluctuation-Dissipation relation 
(FDT) R{t — t') = l39{t — t')dt'C{t — t'). The saddle-point 
evaluation for the two-time variables is the mathematical 
justification of the above-mentioned fact that O, R, D do 
not fluctuate. These saddle-point equations give Eq. (5) 
and (3) (cf. Appendix). 

In this paper we are treating an equilibrium situation. 
Within the liquid phase, this may be achieved by start¬ 
ing from any configuration in the distant past. A more 
practical way, however, is to assume equilibrium at a con¬ 
venient time to (e.g. to =0). How does one deal with a 
non-Markovian equation of motion like (4)? The answer 
is simple: either one makes the memory kernel extend to 
the remote past, or, alternatively, one may assume equi¬ 
librium at to, in other words summing all the past histo¬ 
ries passing through y{to) at to. It turns out [31] that this 
is implemented simply by cutting the memory at a lower 
limit to, as in Eq. (4) (cf. Appendix). This completes 
the derivation of our basic dynamical equations. 


where R)p[x, d] = D[a;, x\p[x, x] and f[xi — X 2 ,Xi — £ 2 ] = 
^-Wlxi-x 2 ,xi-x 2 ] _ I physical p[x,x] is determined 
by 5S/5p[x,x] = 0 and the normalization f T>p[x,x\ = 1. 
The first term in Eq. (13) is an ideal gas contribution 
and the second accounts for interactions. 

Following the thermodynamic treatment [30], we may 
now argue that due to rotational invariance on the hyper¬ 
sphere, p[x{t),xlt)] = p[C{t,t'),R{t,t'),D{t,t')] where 
R{t,t') = {2d/a^)x{t) ■ x{t') and D{t,t') = {2d/a^)x{t) ■ 
x{t'). We can thus make a change of variables in 
the functional integration over x{t),x{t) to Q{t,t') = 
{C{t,t'), R{t,t'), D{t,t')}. The change of variables gives 
for density averages (cf. Appendix): 

(14) 

where strinQ jg Jacobian of the transformation (cf. 
Appendix) and we defined p{Q) = The appear¬ 

ance of the dimension in the exponent leads to a nar¬ 
rowing of fluctuations of correlations, and saddle-point 
evaluation becomes exact (cf. Appendix). In this way 
we can compute the ideal gas term in Eq. (13). For 
the interaction term, that involves two p functions, we 
need the variables corresponding to (xi,xi), {x 2 ,X 2 ) and 


Dynamic transition and timescale separation - We 
now apply the standard MCT methodology to locate the 
density or temperature at which a dynamic transition 
occurs with freezing in a cage, corresponding to the de¬ 
velopment of a plateau in A{t) [1-3, 22]. 

Consider the case when M{t) falls from M(0) to a 
plateau value Mea, and then, at much larger times, to 
zero. Concomitantly, A{t) grows to a plateau value Aeai 
and then continues to grow at a slower (diffusive) pace. 
Denote the fast part SM{t) = M{t) — Mea- In the 
limit in which the plateau times are much larger than 
the microscopic times, and much smaller than the final 
relaxation times, the noise breaks into a fast variable 
5f{f) and a slow random variable (/, as does the friction 
term. Their sum acts as an adiabatically slow field Y(t) 
at those times. We may thus split the equilibration in 
two steps [32]: P (yjF) and Psiow(l^)- When t — t' is in 
the plateau region, we may write the expectations: 


{A{t)B{t')) = J dyP,iow(J") {A)y {B)y , (15) 
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where {•)y = f dF • P(ylV), Peqiy) oc e and 

Peq(y)e = 


P{y\Y) = 


Psio^y) = 


, T -.aP ’ 

/dy'Peq(y')e = y > ( 16 ) 

/ dyPeq(y)e " I “ea j 
•v/27rMEA 


Obviously, Peq(y) = f dV P (yjV) PsiowiV)- We there¬ 
fore obtain the self-consistent equation for Mea: 

Mea = I y dyPsiow(i")(i^)y = X(Mea) ■ (17) 

From Eq. (3) we obtain 


P'^Mea = w— 

^EA 



1 

Aea 


(18) 


The dynamical transition point, at which the plateau be¬ 
comes infinite, happens when Eq. (17) first has a non¬ 
zero solution for Mea = T^/Aea- This point happens 
as a bifurcation, and may be quickly obtained by solving 
Eq. (17) together with M'{Mea) = 0 (cf. Appendix). 

For hard spheres, the result is = 4.807. This re¬ 
sult is fully consistent with the one based on thermody¬ 
namics [12], and in fact Eq. (17) is exactly identical to 
the one that can be derived using the replica method 
(cf. Appendix), consistently with the general RFOT 
scenario [14, 33-38]. This is a particular instance of 
a general correspondence between thermodynamic and 
dynamic results that is verified by the infinite-d solu¬ 
tion, and can be extended to critical MCT exponents [39] 
and to correlation functions [40-42]. In fact, expanding 
around Aea [1“3], one can show that Aea ~ A(t) ~ 
upon approaching the plateau, while A(t) — Aea ~ 
upon leaving the plateau. The exponents a, b satisfy the 
famous relation [1-3, 39, 41]: 


r(l-a)2 _ T{l + b)^ 

r(l-2a) “ F(1 -^26) 


(19) 


For hard spheres, we obtain A = 0.707 which implies 
a = 0.324 and b = 0.629 [43]. Finally, from Eq. (10) one 
can show that the factorization property of MCT [1, 2] 
still holds in d —>• oo, namely that close to the plateau, 

- '(’g.EA-^'/> 9 .EA(A(i) “ Aea) factorizes in a 

function of q and a function of t (cf. Appendix). 

Diffusion, viscosity, Stokes-Einstein relation - At 
long times, in the liquid phase (f < the motion is dif¬ 
fusive and A(t) ~ {2d^D/a^)t, where D is the diffusion 
coefficient. Plugging this form in Eq. (3), and recalling 
that M(t) decays to zero over a finite time, we obtain an 
exact result for D: 

2df T 

0-2 7 -I- /3 dt M{t) ^ ^ 

At low density M{t) = 0 and we recover the bare diffu¬ 
sion coefficient D = T/y. Upon increasing density, M{f) 


increases and D decreases. For ip ^ (p^, where a finite 
plateau of M{f) emerges, J^dtM{t) diverges and the 
diffusion coefficient vanishes as H ~ (^d — ip)'^ with the 
exponent 7 = l/(2a) -I- 1/(25) = 2.34, which is consistent 
with the numerical results of [44]. 

Within linear response theory, the shear viscosity rjs 
of the liquid is given by [19, 45] 

o noo 

VS = dt (cr^i.(t)cr^,.(0)) , 

( 21 ) 

cr^At) = - xpV'^U(x, - Xj) , 

i<j 


where p v are two arbitrary components of the 
stress tensor In the Appendix we show that 

{<yfiu{t)(rnu{d)) = dN M{t), and thus 

^00 

Vs = Ppd / dtM {t) . (22) 

Jo 


This relation shows that for ip —)■ pff , vs ^ IfD ^ (^d ~ 
p)~'^, as it is found in MCT. 

Combining Eqs. (20) and (22) we obtain a relation sim¬ 
ilar to the Stokes-Einstein relation (SER) 


D = 


T 


2d 

pa 


2 Vs 


Tpa^ 1 
2 d T]s 


(23) 


where the second expression holds close to pd- This rela¬ 
tion is interesting because it predicts that the SER is not 
exactly satisfied in dense liquids: the quantity Drjs/T, 
which is constant in SER, has instead a small variation 
proportional to p. This is in agreement with results of [46, 
Fig.7b] that show a linear variation of Dr/g with p in the 
dense regime for large enough dimension. 

Conclusions - In this work we obtained an exact solu¬ 
tion of the dynamics of liquids in the limit of infinite spa¬ 
tial dimension. The picture that emerges has a relevant 
“caging” lengthscale that is smaller than the particle ra¬ 
dius by 0(l/d) (it is about 1/5 for real colloids [47]). The 
physics of diffusion stems from interactions at that small 
scale, while all particle motion beyond that scale con¬ 
sists of uncorrelated steps of displacement, and memory 
of what happens at distances ^ I/d is lost. Diffusion co¬ 
efficients and viscosity are thus decided at distances much 
smaller than the particle radius. This strongly suggests 
that the wavevectors q that matter for this transition are 
not the ones associated with the first neighbor distance 
1 /cr, but rather the much larger ones ~ d/cr correspond¬ 
ing to the cage size. Note that at the transition cages 
are correlated over large distances [33, 48-50]. Also, an 
expansion in the number of collisions [7] seems difficult 
to reconcile with our results, because we expect multiple 
collisions within a cage. 

Is the high-dimensional dynamics related to MCT? 
The answer is that the result is not the one obtained 
from the usual procedure for building up a MCT equa¬ 
tion, which for example gives a different scaling of ipd 
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with dimension [15, 16]. Instead, the system we obtain is 
formally quite close to the slightly more general case of 
soft-spin mean-field dynamics, Eq. (9), because we have 
mapped the system into a one-dimensional dynamics in 
the presence of a colored noise and friction, that have to 
be determined self-consistently through Eqs. (4) and (5). 

In the dense glassy regime we obtain predictions for 
the scaling of the cage radius, of the dynamical tran¬ 
sition density and of the parameter A, that differ 
from the ones of usual MCT [15, 16, 18]. Our results are 
fully consistent with those obtained from the thermody¬ 
namic approach [12, 13, 43], which proves the exactness 
of the REOT scenario [33-38] for statics and dynamics 
in d —)■ oo, as conjectured in [14]. Our results are also in 
agreement with numerical simulations of hard spheres in 
large spatial dimension [44, 46]. 

Interestingly, we find that an ergodic liquid phase of 
hard spheres exists for densities ^ = 4.807. This 

implies that hard sphere packings exist (at least) up to 
^d, and they can be constructed easily through a suffi¬ 
ciently slow compression of the liquid [51, 52]. Note that 
the value of ^d is larger than the best known lower bound 
for the existence of sphere packings, ^ ^ 6/e [53], and 
that it took 20 years to improve the previous best lower 
bound (p'^ 2 [54] by a small factor 3/e. Our calculation 
is simple enough that there is hope to transform it into 
a rigorous proof, along the lines of [55]. This would re¬ 
sult in an improved constructive lower bound for sphere 
packings. 

Euture extensions of this work include the investigation 
of the effect of dissipation [56-59], the study of out-of- 
equilibrium aging dynamics [60], and the study of non- 
perturbative processes in l/d through an instantonic ex¬ 
pansion. The thermodynamic partition function of quan¬ 
tum systems is formally very similar to Eq. (27) and could 
also be studied along these lines. 
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I. FORMULATION OF THE DYNAMICS 
A. Definition of the modei 


We consider an assembly of N spheres of mass m interacting through a repulsive finite-ranged radial pair potential 
V{r). The interaction Hamiltonian is then H = — Xj) with Xi{t) the positions of the particles, a is the 

diameter of the spheres, we note V{fi) = V{a{l + fi/d)) and the rescaled force F{fi) = In the following, we 

define the usual inverse temperature /3 = l/ksT and take ks = 1- In order to exploit efficiently simplifications given 
by the translational and rotational symmetries of the system, we constrain each particle to live on the surface of a 
d+1 dimensional hypersphere of radius R. A particle is thus represented by a point Xi € with xj = R^. 

The volume of this space is V = I]d+i(A) = ^d+iR'^- S and H are respectively the surface and the solid angle. We 
also denote by 


Vd{R) = 


^dl2 

T{\ + d/2) 


R‘^ 


(24) 


the volume of a d-dimensional ball of radius R {i.e. the volume bounded by is Vd+i{R))- Adding this extra 

dimension, rotation and translation invariances in usual d dimensional space are transposed to rotational invariance 
only on the sphere S‘^{R). We recover in the limit R —>■ oo the original definition in a d dimensional periodic cubic 
volume. However, note that in subsections IB and IC, we consider the original model in d dimensional Euclidean 
space, so as to introduce the spherical model only when needed. 


B. The dynamical action 


A possible choice of the dynamics of the particles in the d-dimensional space is the following Langevin process: 

mXi+-fXi =-'^VVixi-Xj)+^i (25) 

^i{t) being a centered white Gaussian noise with (t)Cj' {t')) = 2'yT5ij5^,j5{t — t'). In the following, we will consider 
the overdamped^ case where m = 0. Given a set of initial conditions, the Martin-Siggia-Rose-De Dominicis-Janssen 
(MSRDDJ) [20, 22] generating path integral reads, with Ito convention: 


N 

= / (26) 

N 1,N 

W[xi,Xi,Xj,Xj] 

i<j (27) 

W[x, x,y,y]=i j dt {x - y) ■ VU {x - y) = W[x - y, x - y] 

The sum in Zn is done over all possible trajectories, with the measure given by Dxi = when discretizing 

Gir) 2 

the trajectory Xi(t) in M time steps. Time integrals are taken over an interval noted [tp,fi] in section IV where the 
{a:i(tp)} are fixed (initial conditions, which would correspond to n = 0 in the latter discretization) and the {xi(ti)} 
are summed over (which would correspond to n = M). In section IV we will need averages of the type {xi{to)xi{ti)j 
with to G [fp,fi], and it is not needed to consider times larger than ti in the action owing to causality, as they will 
give no contribution to the average by probability conservation {i.e. for the same reason as Zn = 1). 

The action A is rotation invariant (for both position and response fields using the same global rotation) and translation 
invariant along positions^. Because of the so-called kinetic term, it is not translation invariant along the response 
fields (which are already centered at the origin owing to the white noise), though the interaction term is. 




where the action is: 


A[{xi,Xi)] 

*lx.il = 7/d((rx^ + a.i), 


1 Either Newtonian or Brownian dynamics can be treated by 2 The invariances for response fields follow if one studies the system 
changing the kernel associated to $ in the following sections. stations and translations of the entire system. 




FIG. 1. (a) The different dynamical regimes, described by equation (107): most of the dynamics is determined at the cage size 
scaling 1/d. We assume that the diffusive regime found already at this scale extends trivially at longer times. Concerning the 
slope at very short times, the mean-square displacement is oc if inertia is not neglected; for purely overdamped motions it 
is oc t. Finally, the particles feel the ‘box’ for distances of order R, hence a saturation at this scale, (b) The dynamic triangle 
diagram with three trajectories. 


C. Derivation of the generating functional using a virial expansion 


1. Dynamic virial 


We define S = ™ order to use the Mayer expansion^ as in liquid theory [19]. This ‘grand canonical’ 

form is handy as a generating functional, but note that we assume there is no exchange of particles with a reservoir ; 
InS and IhZn are related by a Legendre transform and virial expansions are more fruitful with the former. We have: 

-too ^ AC „ 

X! + f[xi,Xi,Xj,Xj\) (28) 

i—1 i<j 


with a generalized fugacity z = e * and a Mayer function / = e ^ — 1. We Legendre transform InS with respect to 
iVlnz since one has 


dlnS 

(5(Mnz[a:, x]) 



N 

S{x — Xi)6{x 

i=l 



p[x,x] 


(29) 


where the mean is generated by the functional S. Next, the usual Mayer expansion can be carried out in this 
dynamical case, and inverting the Legendre transform, In S can be written as an ideal gas contribution and the sum 
of all connected 1-irreducible Mayer diagrams"^ with Np[x,x] nodes and f[x — y,x — y] bonds [19]: 


ln2 = — / T)[x,x] p[x^x\{^[x,x\-\-hip[x,x]) + ^ 


A 




(30) 


In infinite dimension, the Mayer expansion reduces to its first term. However, this is strictly true if we assume that 
we are in a regime where the trajectories have the time to wander away only a finite fraction of the box. Because 
we expect (and confirm) that all interesting dynamics (namely, the (3 relaxation with the formation of a plateau 


^ This is why W is cast in a symmetric form so that links in the 
product Ylicj directed. 


4 


1-irreducible means that they do not disconnect upon removal of 
a node. 
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in correlations and the onset of a relaxation towards equilibrium - a relaxation -) occur on such scales, where the 
fluctuation around the initial position is of amplitude 0{l/d) (a scaling consistent with the statics [12, 13, 30, 43]). 
We will show that one even gets the diffusive behaviour which is already decided at the l/d scale (see figure 1(a)). 
To justify this truncation, let us consider for example the second term of this expansion (triangle diagram): 

D[a;, x]B[y, y]D[ 2 ;, z] p[x, x]p[y, y]p[z, z]f[x -y,x- y]f[y - z,y - z]f[z -x,z-x] (31) 


For a finite-support potential V, W{x — y) oc 0(cr — \x — y\) (Heaviside’s step function®), hence 



f[x-y,x-y] 


0 if Vt, 

t[x -y,x-y] if 3t, 


x{t) - y{t) > o- 

- y{t) < ^ 


[0(|a; - 2/1 - cr) - 1]t[x -y,x-y\ 


(32) 


Essentially, t is just a numerical value depending on the details of V (with jrj ^ 2) and the physical content of / lies 
in its finite support, i.e. 


M 

f[x -y,x-y]cz 0(|x - 2 /|-cr)-l = -l-fJ| 0(|a;” -y‘^\-a) (33) 

n—1 


For finite times, a trajectory stays in a bounded region of space, represented as a ball of diameter the typical size 
(X u/d of the trajectory in figure 1(b). The Mayer functions require that each couple of trajectories get closer than a 
at some time. To each set of three trajectories one can associate a corresponding static diagram with three overlapping 
balls (figure 1(b)). One can see this static diagram as the equivalence class of all dynamic diagrams with trajectories 
contained inside these balls. Actually there are lots of trajectories contained by these bounding balls that do not 
contribute because they do not get close enough. Then the sum over trajectories of the value of the integrand is at 
most of the same order of the static diagram value due to the normalization f D[x,x] p[x,x] = 1 which accounts for 
the huge number of equivalent dynamic diagrams compared to the static one. The same program applies to all terms 
in the expansion. We conclude, from [6], that the first term dominates the series® in the d —)■ oo limit. A more careful 
but maybe less intuitive argument is given in IC 2. Note that we did note specialize the discussion to hard spheres 
but to any finite-support interaction potential. The result is then^ 


f TV f 

S = -^ = - V)[x,x]p[x,x]mx,x\+\np[x,x]) +— / B[x,x]Y3[y,y]p[x,x]p[y,y]f[x-y,x-y] 


(34) 


with 


5SI5p[x, d] = 0 from the Legendre transform and the normalization 


additive constants irrelevant for the dynamics. 




. We neglected purely 


2. The Mayer expansion in infinite dimension 

We discuss here the truncation of the Mayer expansion when d —)■ oo in a more pedestrian way. As in IC 1, we 
focus on the triangle diagram in (31) with the Mayer function given by (33). 

Consider two typical trajectories x{t) and y(t) which we expect from the static case to dominate the dynamics at 
large d, and let us focus on the positions x^ and y^ (any other couple would do). These typical trajectories scale 
as follows: the other positions of the trajectory {x‘^}n^i (respectively {//"jn^i) fluctuate around x^ (respectively 
y^) over a neighborhood of size 1/d, see figures l(a),(b). Besides, the typical distance between the two trajectories 
is the diameter a (plus 0(l/d) fluctuations). Then essentially — 2/"| — cr) ~ 0(|a:^ ~ 2/^1 ~ o’) except in 

the neighborhood \x^ — y^\ & [(t( 1 — (j),o-(l -f ^)] where these fluctuations of order 1/d matter {p is of order 1). 


® In the main text it is noted 9 but this notation will be used for 
Grassmann variables from subsection IE on. 

® The critical packing fraction later found in (147) is of order 
0{d/2^) where it is known that the virial series does not converge 
anymore; nevertheless, Frisch Sz. Percus [8] have shown that in 
high d the truncation used here is valid largely above this bound. 


^ In the thermodynamic limit \n3/N = InZjfij/N (as if formally 
the chemical potential was zero in the static partition functions 
analogy). That is why for simplicity we used InZ^^/N in equa¬ 
tion (13) of the main text. 
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However this neighborhood has to be accounted for since its volume V(i(cr)[(l + — (1 — ~ Vd(cr)(e^ —e ^), 

d —^oo 

is comparable to the volume Vd(cr)e“'' of |a;^ — y^ \ & [0, a{l — ^)] where the contribution of f[x — y,x — y] (33) is also 
non zero (for |a:^ — | > ct(1 + ^), f[x — y, x — y] is zero). 

Hence f[x — y,x — y] ~ f*ix^ — y^) with f*{x) being -1 for |x| € [0, (t( 1 — ^)], of order 1 for |x| € [cr(l — ^), (t(1 + ^)] 
and 0 further. What counts for the truncation done in [6] is that '^’d(o') with a around® <7 and /* is 

zero elsewhere. This is the case here since 


/ r 


d—yoo 


-Vd{(j)e ^ 


f 

J[a(i-S).a(i+§)]d 


ocVd{<y){G^^ — e~^^) 


(35) 


Thus we will replace f[x — y,x — y] ~ /Hs(a;^ — y^) = —1 + 0(|a:^ — y^ \ — o'), where /hs is the usual static Mayer 
function of hard spheres, since it does not change the scalings given by f* up to an irrelevant proportional constant. 
The term in equation (31) reduces to 

J D[x,x]D[y,y]D[z, z] p[x,x]p[y,y]p[z, z]fm{x^ - y^)fmiy^ - z^)fmiz^ - x^) (36) 


We can integrate on all variables except x^, y^ and 2 ^ because 

[ ^ j Dai p{x^, ..., x^, x) = p{x^) = constant = (37) 

J (27r)f (27r)i ^ ^ y ^ ^ 

since p{x^) — Sti ~ x))'^ must be a translation invariant function of x^ i.e. a constant. V is the volume 
of the liquid. So the second term in the expansion is actually approximated for large dimension by 

^ J <ixdydzfm{x-y)fmiy-z)fRs{z-x) (38) 

where we defined the average particle density p = N/V. We recognize the same term as in the usual static Mayer 
expansion of hard spheres and can conclude as in IC 1 from [6]. 



D. Spherical setup 


From now on, we constrain the particles to live on the surface of a sphere of radius R embedded in d +1 dimensional 
space, S‘^(i?), see section I A. The field x{t) is promoted to a; : M —and Vt, x{t)‘^ = must be verified. 
Concerning the response field x{t), we will rather constrain it to be orthogonal to the position Held: Vt, x{t) ■x{t) = 0, 
thus living at each time in the hyperplane tangential to the sphere at x(t), cf. figure 2. This way we ensure that 
we will recover rotation and translation invariances for position fields and only rotation invariance for response fields 
once i? —)■ oo is taken^. 

There are many possible ways to enforce these constraints, which are eventually equivalent: 

1. One strategy is to use Lagrange multipliers via a term —Vi{t)xi{t) in the right-hand side of (25), promoted to 
d-b 1 dimensions, which ensures that the trajectory does not get out of the sphere due to interactions or thermal 
noise [22]. It is the one adopted in the main text for concision. 

In d —)■ oo, the value oi Vi ^ -u is non-fluctuating and is obtained by discretizing (25) as follows (in the Ito sense): 


Xi{t -b dt) 


Xi{t) 


-Vi{t)xi{t)dt - -\7xiHdt+ -r]i , 
7 7 7 


iVtVj) = 27 Td%dpi, 


(39) 


More generally bounded by a finite and independent of d quantity 
(for example here it is bounded by 2a for large enough d). 


^ Another way to see it is that the spherical constraint implies 
Xi'Xi = 0 and only tangential fields Xi are needed to exponentiate 
the Langevin equation in the MSRDDJ path integral. 
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At order dt, using that for large d one has B ■ rji —>■ 0 (for any vector B uncorrelated with rji in the ltd sense) and 
rji ■ rji ^ 2dTjdt due to the central limit theorem, we impose the constraint 


= Xt{t + dt) • Xi{t + dt) = Xi{t) ■ Xiit) - —a;i(t) • \vi{t)xi{t) + VXiH] + 

7 7 


and therefore 


/X 1 ^ dT 

+ 


i?2 ' i?2 

We have a general relation [19, Eq.(2.2.10)] for the reduced pressure p: 




For d —>■ oo the fluctuations vanish because we average over d dimensions, we thus have 

^ ^ a:, • ~ ^ x, • = dT(l - p) 

and plugging this in Eq. (41) we obtain that all ^'i(t) are equal and constant in time, given by 


(40) 


(41) 


(42) 


(43) 


Vi{t) 


dT 


(44) 


From the static entropy in d —> oo [6, 8, 12, 29], p = 1 + d<p/2 in the liquid phase, where ip is the rescaled packing 
fraction (p = 2‘^p/d = pVd{cr)/d and the number density p = N/V. 

This choice results in an additional term f dt lyx ■ x in the definition of $ (27), and the summation of paths in (34) is 
over the sphere S'^(i?) for positions and the tangential hyperplane. In practice, we will integrate on the whole d + 1 
dimensional space and enforce the constraint through Dirac deltas. Note that, with x (resp. x) representing the 
position (resp. response) field at some time and E = Span(a;), 


f dx d{x^ — R^) = rtd+i [ dr 

./B'i+i Jo 


R'^ 


V 


1 


r‘^6(r^ — R^) = fld+i— = — = — / dx 
\ ) <^+^2R 2R 2R a. 


f dxS{2x-x)= I dx//dx±5{2\x\x//) = f dx//dx± S{2Rx//) =-— j da;_L 

dR'i+i JexE-^ JexE^ Jex 


(45) 


These choices rescale the path integral measures with respect to the ones on which does not affect the dynamics. 
In the thermodynamic limit of infinite radius R (with p = N/V fixed), we recover the original d-dimensional space. 


2. In the previous method, we would exponentiate the Dirac delta functions, giving additional terms in the exponent 
J dt px-x and J dt jj,{x‘^ — R"^) . We see that the Lagrange multiplier would just shift p. Hence, another technique 
is directly (and somewhat physically blindly) to promote the path integrals in (34) to d -|- 1 dimensions and use 
Dirac deltas to constrain the x, x fields. This is what we are going to follow with fields v, v but note that it is 
not exactly the Lagrange multiplier^ even if it plays a similar role. 

3. We might as well add a soft constraint A'Y/iix/ — R?)^ in the Hamiltonian. It would add a single-particle term 
2A f dtix ■ x(x^ — R/^) to 4>. Then we can write (still at a given time) 

^-2Aix-x{x^-R^) ^ f p^ps(^-X-x) OC f p^J,p^Pe^^’>/A-2^ux■x-iHx^-R^) 

J^d+i \2A J jRd+i 

and in the hard limit A —)■ oo this is the same as enforcing the constraints through S{x'^ — R^), 6{2x ■ x). 
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E. Superspace notation 


1. Translation of the dynamics into superfield language 


As a compact way to write dynamical equations, we will use superspace notation [61, 62], At any step, one can 
unfold this notation to recover the standard dynamical variables. {0i,9i} are Grassmann variables^°. Let us define 
X = ix for convenience. We encode the position and response fields [x^x] in a superfield^^ x{a) = x{t) + 9i9ix{t), 
where arguments are denoted by a = {9i,9i,t). The Mayer function and the kinetic part can be explicitly written: 


f{x) = e-/daE(=r) _ 

®W = 7/d<.^- 


1 with 


d9i 


I da 

iH ) 


ddid^idt 


(47) 


a;(a)^ = x{t)‘^ + 29i9ix{t) ■ x{t) implies for the constraints that 6(x(t)^ — R^)5{2x{t) ■ x(t)) = S(x(a)'^ — R^). 

The measure D[x,x] is replaced by Da; = D[x,x] where integration over x is on the ‘imaginary axis’ Then the 

action can be written in the form: 


S = - J Rtx 5{x{af - R^)p{x){\np{x) + 4>(a;)) + y J BxByS{x{af - R^)6{y{bf - R^) p{x) p{y) J {x - y) 


(48) 


still with 5S/5p{x) = D 


and 


J Dx5{x{a)'^ — R^)p{x) = 1 


2. Definitions about superfields 

Here we provide some useful generalization of usual field-theoretic tools and other identities in superfield language, 
used in the following sections. Proofs are not given here, but are easily obtained using the definitions and by direct 
computations. 

1. Superfields: We will consider (one-component) superfields and define them as x{a) = x{t) + 9x9xx{t). 

We will also use operator (two components) superfields analogous to the replica case, such as q{a, b) = x{a)-x{b). 
Similarly, an operator superfield r can be cast in the canonical expression with Grassmann variables and real 
scalar fields: 

r(a, b) = ri{t,t') + 9i9ifi{t,t') + § 292 X 2 ( 1 , t') + 9i9i9292r2{t,t') (49) 

2. Dirac deltas: For superfields, they are simply functionally defined as 5(x(a)) = 5{x{t))5(x(t)). For operator 
superfields, they are simply functionally defined as a product of the functional deltas of their components 
appearing in notation (49). If the superfield is symmetric we need to introduce deltas only on the independent 
part, which is the case for q. 

3. Path integral measure: We clarify here the path integral measure for future needs. For a general superfield r, the 
path integral measure is defined as Dr = D[ri, fi, f 2 , r 2 ]. For symmetric superfields such as q, we will only sum on 
the symmetric components of it, and call it D^q = D*[(7i, g 2 ]D[( 7 i, 92 ] with D^gi = Dqi Htx' 5(di{i, t') — Qiit', t)) 
in a discretized point of view. Therefore q(a, b) with ‘a > 6’ (loosely speaking) will appear in the path integral, 
the previous Dirac deltas imposing them to be < 7 ( 6 , 0 ), so that we can use all components and introduce q as a 
symmetric superfield in the path integral as it should be from its definition. 


Let us emphasize here that this is only a compact notation with¬ 
out any physical content (for our present purpose), but we note 
by way of excuse that the computation is prohibitively compli¬ 
cated proceeding otherwise. 


We do not need to resort to purely fermionic components (i.e. 
linear in 81 and 81 ) in the superfield x since we interpreted the 
Langevin equation in the Ito sense, implying the MSRDDJ Ja¬ 
cobian is 1. Any other discretization scheme can be equivalently 
considered by adding two such ghosts. 
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4. Product of superfields: We define the product of two superfields r = pq as the generalization of the operator 
product 


r{a,b) 


dcp{a, c)q{c, b) 


(50) 


5. Identity: The identity for superfields reads in components 

l(a, b) = 9i6i5{t — t') + 9292S{t — t') (51) 

6. Inverse: The definition of the inverse r~^ in the superfield sense reads, through the relation r~^r = rr~^ = 1, 

r~^ = /3 + 9i9if3a + 026*27/3 + 0i0i0202(7/?Q; + rf^) (52) 

with fi = {r 2 — rirf^f 2 )~^, 7 = —rf^r 2 and a = —rirf^. 

7. Determinant: We define a ‘superdeterminant’ for superfields as 

sdet(r) = det(ri) det(r 2 — ■firC^f 2 ) (53) 


8. Trace: The trace is defined by str(r) = / dar(a, a) = J dt (fi + r 2 ){t,t). 

9. Integral representation of Dirac deltas for superfields: an integral representation is expressed as 


ita) = / 

6{x) = J 


Dp e* ^ v{a,b)q(b,a) = J Dp e*"tr(pq) for a superfield operator 
Dy similarly for a one-component superfield 


(54) 


For a symmetric superfield such as q{a, b) = x(a) ■ x(b), we only need to introduce the independent part of it, 
that is, taking the additional superfield p to be also symmetric, we only have to sum in the exponential over 
half of / dad5p(a, b)q{a, b). Rescaling p with the factor i, the formula above is unchanged^^ provided that the 
measure is understood as D®p since p is symmetric. 


10. Derivation with respect to a superfield: Defining the derivative with respect to a superfield in a way similar to 
operators, using the convention for functional derivation: 

= (55) 

where r is a superfield, and 1 is the equivalent of the Dirac delta for superfields. Similarly for a symmetric 
superfield q we get 

= Ha,c)l{b,d)+ l{a,d)l{b,c) (56) 

dq(c, d) 

which gives the factors 2. We will thus use the formulas 

= (57) 


11. Gaussian integration on superfields: A direct calculation with components shows that the Gaussian integral can 

f dadb q~^ (a,b)h(a)h(b) 


be cast in the familiar form: 

[ dadb q{a,b)x{a)x{b)-\-J dah(a)x{a) _ 


Csdet(q) 


(58) 


h and x being here scalar fields for the sake of clarity {i.e. d = 1 ). <7 is assumed symmetric, i.e. qi and <72 are 
symmetric operators and qi = qj■ The Gaussian integrations require that qi (respectively ( 72 ) and <72 — qiqf^q 2 
(respectively qi — ( 729 ^^ 91 ) are positive definite. 


12 


All numerical constants, such as this one coming from the rescal¬ 
ing, will be omitted as we will eventually calculate all propor¬ 


tionality constants in another way. 
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12. Gaussian moments: Moments of a Gaussian random variable in notation SUSY enjoy similar properties than 
discrete Gaussian random variables. As an example, we can compute the following variance: 

Da?a;(a)x( 6 )e- 5 /dad 6 q(a.b)x(a)x( 6 ) ^ q-i(a^ 5 ) ( 59 ) 



13. Other useful identities: The following relations are useful for the derivation of the interaction term, for any 
superfields A and B such that the series converge: 


lnsdet(A + B) =lnsdetA + stry - 1 -= InsdetA + str V - 1 - 






(60) 


(A + B)-^ = y (-1)"(A-1B)" A-i = A-i y (-1)”(BA-1) 


k n>0 


n>0 


These formulas are just generalizations of power series to superfields. The radius of convergence of the series 
used is 1 . 

14. The projector P: We define a superfield P{a,b) = 1 and A a scalar. We get the following important formulas: 


(XP + B)-^ = B-^ 

Both formulas require |str(PB' -i)|<l/|A|. 


lnsdet(AP + B) = Insdet(B) + In |^1 + Astr(PS 

A 


1 + Astr(PP 


-1 


PB 


-1 


(61) 


II. TRANSLATIONAL AND ROTATIONAL INVARIANCES 

We now take into account rotational and translational invariances and take the limit d —>■ 00 . In some cases the 
order of the two limits is irrelevant, but when relevant, we should take the i? —>■ 00 limit first. In other words, we 
should consider for example that R/d is a large quantity. 


A. ‘Functional spherical coordinates’: invariances using the mean-square displacement 

As emphasized in section I A, the aim of the introduction of S‘^{R) is to take into account both translation and 
rotation invariances on Euclidean d dimensional space by only rotation invariance on a sphere of a d + 1 dimensional 
space, which is actually easier to handle in the viewpoint of the dynamics. Indeed, the MSRDDJ action A in (27), 
now promoted to d + 1 dimensional fields, and the constraints in (48) are invariant by the same rotation TZ for both 
fields i.e. (x,x) —> (7Zx,7Zx), which is transposed to superfields as a global rotation^^ x —> TZx. 

Now considering expression (48), we define a superfield q(a, b) = x{a) ■ x{b). We assume that the d + 1 dimensional 
liquid is invariant by rotation i.e. p{x) = p({x(a) ■ x{b)}a,b) = piQ)- This way we will remove all irrelevant variables 
and be able to use a saddle point method. 

Eventually, as regards the i? —)■ 00 limit, it is more convenient to consider the mean-square displacement (MSD) 
D{a, b) = {x{a) — x{b)Y since it is a finite quantity as long as the difference t — t' is finite, at equilibrium. Before the 
dynamical transition is met, D is of order when t — t' —>■ 00 . We thus expect an artificial second plateau of order 
R^ due to finite size effects, that will be removed when R —>■ 00 , giving back diffusion at long times (see figure 1(a)). 
One can check explicitly that in this limit the original d-dimensional MSD is recovered and that it is translation and 
rotation invariant for the position and rotation invariant for the response field, see figure 2 . 
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See subsection ID and footnote 2. 
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FIG. 2. Notations on the sphere for x, x. The vector x{t) — x{t') is along a chord and when R —>■ oo, it lives in the same 
d-dimensional space as x. 


B. Scalings in the infinite d limit 


First, from now on we measure distances in units of the diameter ct, i.e. we take a — 

To use the saddle point method, we need to specify how the quantities defined here scale with dimension. From the 
statics in [12], if we define u = x — X where X denotes the translational degree of freedom of a particle so that u 
characterizes motion around an average position, u scales as 1/d. As $ is translation invariant along x. 




(62) 


ft ^ 1/d implies that, for the action to describe the Langevin process above ((25)), we need the two terms in the 
integral to scale identically in the limit d 


1/d as well. This implies that D — A./d 


oo. Thus u 

order 1. Hence, for convenience, noticing that D = 2B?P — 2q, we dehne the rescaled quantities 


where A is of 


Aiiq = 2dR^ 


Q = AiiqP — A = 2dq 


(63) 


Aiiq corresponds to a typical MSD between particles (in the liquid phase) on the sphere S‘^(i?). Q{a,a) = A^q is the 
spherical constraint. 


We then choose 


7 = 2jd^ 


in order to write 


= dd) with 7 and of order 1 so that $(a;) scales like lnp(a;), otherwise 


the infinite d limit would not be well defined. 

Indeed, from the statics of the liquid phase [6, 8, 12, 29, 30], we will assume that p is exponential in d, that is 


p(q) = A(g)e^^(«) i.e. 


Inp(q) ~ dl}{Q) 

d—¥(x> 


(64) 


with Q of order 1 and A a subdominant factor {i.e. non-exponential in d). 

Concerning the potential V{r), r will typically scale here as r = a{l + where p, is of order 1. Besides the two terms 
in the action A in (27) must have the same scalings to describe the underlying Langevin process (25). As A can be 
cast in a form similar to (34) using p[x, x], thus similar to an ideal gas term plus an interaction term with ‘propagator’ 
FF, the following derivation of the infinite d limit will indicate that such interaction terms have the same scalings as 
the ideal gas term only if the ‘propagators’ W (or /) scales as 4) i.e. 1 here. Thus V, which scales like W (cf. (47)), 
must be of order 1. For example, the soft harmonic spheres potential is Vss(r) = k( 1 — r)^0(l — r) = K{p/d)'^Q{—p). 
The amplitude must he k = with k of order 1. 
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C. Ideal gas term 


We write equation (48) as 5 = 5 ig +5int and focus on the ideal gas term 5 ig- Exploiting the invariances, we apply 
the program mentioned in II A: 

5ig oc J r)x S{x{a)^ — R^)p{x){\np{x) + ^{x)) 

= f DxD^^qSiqiujb) - x{a) ■ x{b))6{q{a,a) - R^)p{q){\np{q) + ^{q)) 

% (65) 

= / DxD"qD'^q',5(q(a,a) 

= / D'^qD'^q' 6{q{a, a) - i?2)e*"f(‘?9')-4i insdet(2*q')p(q)(inp(q) + $(q)) 


The explicit expression of ^{q) will be computed in III A. The integral over q' is evaluated through a saddle point 
method for d —)• oo. Using the tools from subsection IE2, the saddle-point equation is 2iq = (d + l)q'~^|sp. We 
introduce the rescaled variable Q and neglect subdominant terms which will be calculated in subsection IIC. For 
convenience with respect to the saddle-point equation in section III, we exponentiate the resulting spherical constraint 
Q{a,a) - Aiiq, 


5iG OC J D[Q, u]e^ Insdetg-f /dai.(a)(Q(a.a)-A„,)^^Q)(jj^^(g) ^ 


( 66 ) 


In the limit d —>■ oo we apply a saddle point method, thanks to the scalings provided in IIB. We can get rid of 
proportionality constants using the normalization of the density: 


__5ig_ 

J Da? 6(x(a)^ — R^)p(x) 


/D[Q,».]C(g)eir(Q.^)(lnp(Q) + 4>(Q)) 
/ D[g,iz]C(Q)e5r(Q.^) 


(67) 


where C(Q) accounts for forgotten subdominant contributions to the integral, and 


r(g, ly) = Insdetg 


da iy(a) (Q(a, a) — Auq) -|- 2fl(Q) 


( 68 ) 


is the saddle point function. In d —> oo we maximize T, in particular 


dr 


sp 


giving the result 


(5r 

= 0 and — = 0 Q(a, a) = A^q 

oiy sp 


5ig = -d[r!(Q^p) + $(g^p)] 

d—¥oo 


This can be expressed explicitly using once again the normalization condition 

J D[g,z^]c(g)e5P(Q’^) = I 


(69) 


(70) 


(71) 


Evaluating it for d —)■ oo and taking the logarithm of the resulting equation at dominant order 0(d) provides 


r(g®P, i/^P) = 0 up to irrelevant additive constants^"*, from which we get 


5ig = ^ Insdetg^P - ^ /daiz^P(a) (g^P(a,a) - Auq) - d$(g^P) 

d—^oo Z Z J 


( 72 ) 


These constants were hidden in C(Q®p) and do not depend upon 
therefore having no influence on the dynamics. 
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D. Interaction term 


1. Changes of variables to ‘funetional spherical coordinates’ 


Now focusing on 5int, we introduce the superfields q(a, b) = x{a) ■ x{b), p{a, b) — y{a) ■ y{b) (through a symmetric 
measure D®) and the interaction superfield uj(a) = (x(a) — y(a))^, a one component variable since the Mayer function 
f[x — y,x — y] in (47) needs only scalar products between its variables at equal times. 


5int oc J D[a;, y](5(a:(a)^ - R^)5{y{bf‘ - R^)p{x)p{y)f{x - y) 


= I D[=r, y, q, q',p,p', a;, a;']<5(q(a, a) - R^)5{p{a, a) - p{q)p{p) f {^) 

^ g-* / dadb [q' (a,b)x{a)-x(b)+p' ia,b)y{a)-y(b)]-i J da lj'{ a)l{x(a)-y{a))^-ui{a)] 

= J D[g,9',P,p',‘^,‘^']'J(9(a,a) - R^)d{p{a,a) - R^)p{q)p{p)f{^/^) 


X exp 


^istr(<7q' + pp') + i 


da u)'(a)u:(a) 


d + 1 
2 


In det 




(73) 


The last line is obtained by Gaussian integration on superfields x and y, introducing the superfield 

u){a,b) = u)'{a)l{a,b), and det means the determinant of the 2x2 block matrix consisting of superhelds. We now 

change variables to exploit the x ^ y symmetry in (73), 




Q±P 

2 


and 


g± = 


q'±p' 


(74) 


whose Jacobian is subdominant. The 2x2 block matrix can be transformed through a (Weyl) rotation which does 
not affect the determinant: 


/l -l\ fq' +u> -u) \ f 1 l\ _ fq'^+20 q'_'\ 

71 1^1 p' + d;j 7ll^-l q'_ q'J 


(75) 


In d —)■ oo, owing to the x ^ y or q ^ p (respectively q' 44 p') symmetry in (73), the saddle-point value of q_ 
(respectively q'_) is zero. Dropping the J- index for the two other superfields^^, we get 


5int OC J D[q,q',u;,u;']J(q(a,a) -i?2)e2*str(qq')+z/da^'(a)u.(a)-4il„sdetq'-i*±il„sdet(q'+2.S)^(g)2^(^) 

To simplify the last term in the exponential, we make the following change of variables: 

q = R^P-^D and q' = (^R^P - (77) 

The Jacobian is subdominant ; all such subdominant terms will be calculated in subsection IID 2. The aim is to 
recover the same saddle point function T at 0{d) in the exponential in 5int as in the ideal gas term. Let us focus on 
the last term in the exponential, neglecting irrelevant constants: 


D' 


lnsdet(q' J- 2u;) = — Insdet I R?P -— + Insdet 1 — - — -coD' + In I J- -— -str 


2i 


D' 


= — Insdet 1 R P -) J- Insdet | 1-,- coD ) J- In str 

a J- 1 


d+1 

2i 


2R^ 


d J-1 


D' 


2i d J-1 


D' 

2i d J- 1 


(78) 


The symmetry also implies that the saddle-point value of q_|_ is 
the same than q or equivalently p (same for the primes). 
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We used (61) and expanded using the limit R 


to this section, which is valid if str 


2i d+l ) 


oo before d 
-1 


oo in the last line as emphasized in the introduction 
is not zero. 


We expect the same scalings for D and D', so let us assume D' = A'/d and u>' = dfi', hence usD' = /xA' = 0{d^). 
Using (60), we can use expansions: 




In str 


D' 

‘ 2/1 d “t" 1 


-1 


, str(PuA'/x) ^ / 1 

= lnstr(P/x) + ——- , +0 ( ^ 

d-j-oo d+l str (Pfi) \d^ 


Summarizing, 


(79) 


X [ D[a;, / da In str(J’/x) —'istr(f’/xA^/x) / str(f’/x) 


(80) 


At order 0(d) in the exponential in 5int, we get (twice) the same terms that we had in (65). ; therefore in the last 
line, all dependence in D' is of 0(d°). We now go back to the original variables q and q' by making again the change 
of variables (77). The terms in the exponent that must be kept for the saddle point in d becomes 

2istr{qq') — (d + 1) lnsdet(2iq') + 2 In p{q) (81) 

which is exactly twice what we had for 5 ig. Hence, 2iq — {d+ l)< 7 '|gp~^ i.e. P'bp = D. 

Then the term str(/2A) = J da fi'(a) A {a, a) = 0 because of the constraint. We also note that 

u;(a) = {x{a) - y{a)f = {x{t) - y{t)f + 2di6»i(a;(t) - y{t)) • {x{t) - y{t)) (82) 

As we expect {x{t) — y{t)Y = ct^( 1 + 0(l/d)), we define u}{a) = 1 + 2/x(a)/d. We set A = / day,'{a) = str(P/i) 
with a Dirac delta and exponentiate it with a conjugated as usual. The interaction term^® now reads, once again 
exponentiating the constraint S(Q(a,a) — A^q) through a superfield iz. 


‘5int oc J D[Q, iz]e‘^^^‘^’‘^^P(( 5 ), with P defined by: 

P(Q) oc J D[/x, /x']dAdA' exp ^zAA' + idX — ^ ^ ^ In A + 2 z J da/x'(a)(/x(a) — A'/ 2 ) 

X exp (—zAAiiq + istr{Pp,Qp,)/X) / 



(83) 


We used the simplification str(P/xP/x) = A^. The integral over A can be performed in the infinite d limit: the 
saddle-point equation gives A|sp = l/2i. Performing the following steps: 

1. In the Mayer function /, expand the square root yd^F2/x(a)7d in the limit d ^ oo 

2. Rescale y{a) — A'/2 —s- y{a), still calling the new superfield y 

3. Rescale A'/2 —> A, dropping the prime for convenience 

4. Perform the Gaussian integration on y' 
we obtain 


T{Q) oc 


e-Aiiq/2 

-y/sdetQ 



+ 


fJ'd- X \ 

d J 


(84) 


For now we put all subexponential dependence in + ; explicit 
expressions will be given in IID 2. 
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2. Normalization 


Note that all the nonexponential in d dependences overlooked during the procedures of the different changes of 
variables does not depend upon the choice of the Mayer function /. Here we benefit from this to give the explicit 
expression of 5int- 

In the MSRDDJ action A (27) we sum on times belonging to an interval [tp,ti], where initial conditions are fixed at 
tp. ti labels the final state, and if we sum on all positions at ti, we have Zn = 1 in (26). Let us pick s G]tp,ti[ and 
define a test function fo[x,x] = ©(1 ~ 13^(5)|)- Note that the choice of the test function is not completely arbitrary: 
as seen in IC, it should satisfy the properties of the true Mayer function / that we used to derive S, as it must reject 
all trajectories that do not get close at some time. Making a choice that does not respect these properties would lead 
to absurd results. We obtain, using first the expression of 5int in (34) and setting y = u + x and y = u: 


‘5int[/o] 


N 

Y 



, x]D[u, u\p{x, x]p[u + X, u]0(l 


u{s) ) 


N 

~2 


N ( (27r)5 Y V 

2 V j (27r)f 


du"= 

(27r)5 


0(1-K=|) 


pMi) 

2 


dx"- du"° 

(27r)5 {2TT)i 


p{x-^)p{x^‘ 


+ u”=)0(l-|w”1) 


(85) 


As in IC 2, we discretized the trajectories and used that translation invariance and the normalization f D[x, x]p[x^ x] = 
1 imply / Dx d p[{x^i = constant = (27r)5/V. Us labels the time s, i.e. s = tp + ns{ti —tp)/M. 

(27r) 2 

We define C'{Q, A^q) accounting for all the overlooked terms. We have, taking the saddle point over Q, 


Ynt[/] ^int[/] 

Sint[fo] pVd{l)/2 

e-A>io/2/dAD[(5,iy,/x]C'(Q, Aiiq)e''r(Q.‘^)e^-5/dad6M(a)Q-ba.b)M(fc)j (l + 
e-Aiiq/2 J dAD[g,IZ,/x]C'(Q, Aiiq)e'^r(Q,i.)gA-i/dad6M(a)Q-i(a.b)M(fc)jjj 

= _ \ _ f p)udAe^“5/dod6M(a)Q“^bp(o.&)M(b) f 

C"(Q‘^P,Aliq) 7 ^ \ d J 


C" is given by 


C"(g"P, Aiiq) = y '|sp(».f>)M(fc)g) 

= J _ gAiiq/2 


( 86 ) 


(87) 


where we introduced the superfield g{a) = —6i6iS{t — s) to integrate the Gaussian. We also used the constraint 
Q’'P(a, a) = Aiiq. Now we can conclude: 


where J'iQ) 


Sint — 


pVd(l) 


d — 


j-(g^p) 


e-'^iiq/2 

•y/sdetg 





( 88 ) 


E. Final result in the limit d —> cxa 


Collecting the results from the last two subsections, and using (61), we obtain the final result in the infinite 
dimension limit: 


S = Ynsdetg-^ [ dai>{a) (Q{a,a) - Anq) - d^{Q) + ^iF{Q) 
d^oo 2 2 J 2 


sp 


(89) 


up to irrelevant additive constants (cf. IIC), with if defined in (88). 
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III. SADDLE-POINT EQUATION 

A. Explicit form of $(Q) 

We now make explicit the Q dependence of <I> using <I) justified in subsection IIB. From (27), it is Gaussian in x 
and X, 


^{x) = ^^{x) =jd Jdt^x-x — Tx^^ =2d J dtdt' x{t)k{t,t')x{t') + 2 j dtdt' x{t)k{t,t')x{t') 

f) 

with k{t, t') = —^TS{t — t') and k{t, t') = ^ ~ t') 


(90) 


The kernel k is symmetric while k is antisymmetric. Let us define a symmetric superfield: 

fc(a, b) = k{t, t') — 6i6ik{t, t') + 6292k{t, t') 


(91) 


One can check that $(<5) = str(feQ) 


gives back expression (90). 


B. Saddle-point equation for the dynamic correlations 


From subsection ID, the probability density of trajectories p is given by the saddle-point equation SS/Sp{x) = 
SS/Sp{Q) = 0. In the infinite d limit, S depends on p, or equivalently on its logarithm Q, only through its saddle- 
point value 11((5®P). From the relation F((5 ®p,i/®p) = 0 derived in IIC thanks to the normalization of p, Q((3®p) is 
explicitly determined by the saddle-point values Q®p and iz®p. Hence the saddle point condition is equivalents^ to 
(55/(5Q|gp = 0 in (89), and is, as a consequence, equivalent to the saddle point condition used in the virial terms for 
d —)■ oo. 


6S 

SQ 


sp 




_^l + 2fc + |^ =0 

2 oQ sp 


(92) 


The derivative of T is 

= J^{A)Q-'^{a,b) - J da'db' Q~^{a,a') J J Vfj, fi{a')n{b')f (^1 + 


/X + A 
d 


Q-\b\b) 


with the Gaussian measure / T>fi • = 


1 


•y/sdetQ 


D/x • e 


— dadb fx(a)Q ^(a,b)fj.(b) 


(93) 


Hence the saddle-point equation, V(a, b), 



(94) 


Together with the spherical constraint Q^p = A^q which shall provide iz^p, this determines Q^p. 


C. Simplification of the saddle-point equation 

1. Exploiting Ward-Takahashi-like identities 

Here we drop the labels ’sp’ for convenience. Generically, derivatives of S are needed for example to compute the 
Mode-Coupling Theory (MCT) exponents. They can be simplified using Ward-Takahashi-like identities. From the 


The condition (5iS/<5t'|sp = 0 is once again the spherical con¬ 
straint Q®P(a,a) = Aiiq. 
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definition of the Mayer function / (47), a quantity like (•/)^ is a difference between two averages, one with potential 
V and the other without. Let us focus on the non Gaussian part, with potential: 

J = (1)^ (95) 

Usually one shifts /x by e in some well chosen average {0)y and demands the cancellation of all non-zero orders in e, 
but a more compact way to get moments of Q/jl is to write, V(a, &), 


-y/sdetQ J 

= Q"i(a, b) + (F(/x(a) + A)U(/x(6) + A))^ + (U'(/x(a) + A))^ l(a, b) 


where we integrated by parts twice. This method can be easily generalized to higher moments. 


Q I (a,6) 




6fj,{a)Sfj,{b) 

- /dcy(p.(c)+A) 


+ Q ^{a,b) 


o—5 / dcde/i(c)Q ^{c,e)fi{e) 


(96) 


2. The value of T{Q) at the saddle point 


The measure in (95) can be interpreted as an average over a Langevin process with potential V. Provided we sum 
over all possible trajectories of fi, equation (95) is actually the normalization of probability (l)y = 1. 

Similarly fF{Q) can be interpreted as a difference between averages over two dynamical processes, one with potential 
V and the other free. 


J-(Q)= J dAe^-'^'^'/^y = J (j;^gA-Au,/2 J 


VpLe 


— f daV (/x(a)+A) _ 


Vfil 


= J dAe ^-^'-/2 [( 1 )^ - ( 1 ),] 4 » | j -( Q ) = 0 


(97) 


Hence T is zero for all acceptable dynamical propagators (positive definite, as we expect Q is, at least at its saddle- 
point value dominating the dynamics) due to normalization. 


3. Definition of the memory kernel M 

From equations (94), (96) and (97), the saddle-point equation is simplified as: 


Q ^{a,b) = 2k{a,b) — M{a,b) -b (i>'(a) + diA(a))l(a, &) 


where 


M{a, b) = ^JdX + A)U(/x(6) + A))^ 

= / dAe^-^”-/2^U'(/x(a) + A))^ 


(98) 


(99) 


In our study of the dynamics of the system, M will play the role of the analog of the MCT kernel. 


IV. EQUATION FOR THE DYNAMIC CORRELATIONS 
A. Equilibrium hypothesis 

In this subsection and the next ones we will unfold the SUSY notation to get rid of it, coming back to the standard 
dynamical variables. 
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We focus on the equilibrium dynamics of the system, assuming that time translation invariance (TTI) as well as 
causality hold, and consequently fluctuation dissipation theorem (FDT): we assume we start in the remote past (time 
tp, formally sent to —oo), so that the system is at equilibrium when a finite tg is reached. 

Q = 2dq (equivalently A) takes its equilibrium form and so does M: 

Q{a, b) = C(t — t') + 9i9iR(t' — t) + 9292R{t — t') 

( 100 ) 

Mia, b) = Mit - t') + 9i9iMit' - t) + 6 » 26 » 2 M(f - t') 

Along with TTI and (98), this implies that i>'(a) = i/ and di>'(a) = Si' are constant and real quantities. From this the 
inverse Q~^ reads: 

Q-\a,b) = d{t-t')+9i9iRit'-t) + 9292R{t-t') with C = -R-'^CR-^ and R = R-^ (101) 

C and R satisfying fluctuation-dissipation theorem, one can check that C and R also do. Similarly, we can verify 
directly with (90) that k and 1 verify the latter relation. We conclude, from (98), that M also satishes FDT (which 
was suggested by (99)). 


B. Mode-coupling form of the saddle-point equation and the effective stiffness 


We can cast the saddle-point equation into a mode-coupling form by multiplying (98) by Q on the right. The scalar 
component of the equation obtained reads at 

0 = - 2{kC + kR^) - MC - MR^ + {v + 5v)C 

=jdt>Cit' — t) — 2jTR{t — t') -\- {v -\- 6v)C{t — t') — P J drt (9„M(f' — u)) Ciu — t) (102) 

— /3 / duM{t' — u)duCit — u) 

J —OO 

We can assume, for instance, that t > t', and use FDT: 

0 =^Cit — t') + (v + Si')Cit — t') — l3 j du M{t'— u)duC{t — u) — P f dudu[Mit'— u)Cit — u)] (103) 

J t' J —OO 


Using the relaxation for long times and making the substitution v = t 1' — u, 


7 C(f -t') = -{v + Sv- /3M(0)) Cit -t')-P f dv Mit - v)Civ - t') 

Jt' 

t' —)■ t~ gives v + Sv — /3M(0) = 

Aliq 


(104) 


(7(0) = —TRi0~^) represents the immediate response to a perturbation. At such very short times the potential is not 
relevant, particles follow a free dynamics. (7(0) can thus be obtained using (25) yij = with vi = dT/Rp 

since, at equilibrium, TTI and the equipartition theorem holds. One can solve this Ornstein-Uhlenbeck process and 
compute (7(f) = A^qexp(—dTf/yi?^), giving —7(7(0) = T, hence 


v + Sv — /3M(0) = 

Aliq 


We conclude that the mode-coupling-like equation for C is, for t > t': 


jCit — t') =—-^—Cit — t') — P [ dv Mit — v)Civ — t') 

Aliq Jt' 


or equivalently for the MSD A = d D = A^q — C at t > f: 



, T 

/■* 

pA-it — t 

') = T-—Ait-t' 

Aliq 

) — P dv Mit — v)Aiv — t ) 

Jt' 


(105) 


(106) 


( 107 ) 
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C. Effective Langevin process 


The aim is to compute M, the mode-coupling-like kernel, as a function of to solve the saddle-point equation 
for providing correlation and response of the system. To do this, we must calculate correlations of the force F 
at two times to and ti ^ to- To achieve this program, we will interpret, as mentioned in IIIC2, the average defining 
M as two-point correlation functions of a Langevin dynamics with potential V. We drop the notation V in averages 
since we will be refering only to this process from now on. To this end, let us unfold the MSRDDJ path integral in 
SUSY notation (95) using the saddle-point equation (98), where /x(a) = /r(t) -I- 


^ ^ -v/sdetQ 

1 

VsdetQ 


D/Xe “5 f dadb/j.(a)Q (a,b)/j.(b)-f da V(/J.(a)+^) 

D/Xe~^ / dadb fi{a)]^{u-^Si^)l-\-2k — Al'^{a,b)fi(b) — f da y (/i.(a)+A) 


(108) 


VsdetQ plays the role of normalization^®. The corresponding Langevin process with potential V, depending on A, is: 

= + [ dt'+ X)+ C{t) 

■ftp (109) 

with (^((t)')=0 and = 2'^T5{t — t') F M{t — t') 

We used that M is causal and consider times t € [to,ti\. M (and M by FDT) can be computed self-consistently with 
force correlation functions of this process through its definition (99): 


Mit - ^') = I / dAe^-^'-/2 {F{fd{t) + A)F(/t(t') + A)) 


( 110 ) 


Using FDT for M, an integration by part with the fact that t — tp > to — tp ^ Ta the relaxation time of the system, 
above which correlations vanish, and (105), we get the generalized Langevin equation equivalent to (109): 

iKt) = - P [ dt'M{t-t')fi{t')+F{^i{t) + X) + Cit) (111) 

^liq Jtp 

Note that so far in this section we have used equilibrium properties for all observables. It is either because we 
considered them at a time t ^ to where equilibrium is reached, or in the case of convolution products with M where 
the integral extends to the remote past, because we expect that M (respectively M) vanishes quickly (respectively 
vanishes on a hnite time scale Tq,). Then for hnite times where they do not vanish, the system has equilibrated from 
the initial condition in the remote past tp, and we can consider their equilibrium properties. 


D. Memory kernel in equilibrium: resumming trajectories from the remote past 


As emphasized before, dynamical two points functions at (toUi) should be function of a single argument, the time 
difference ti — to, as we assume that equilibrium is reached. But equilibrium properties hold only if we consider times 
t ^ to- Therefore, we wish to use to as our initial time, but in the understanding that at this time the system is at 
equilibrium. Similarly, we are not interested in what happens after U, as it should be irrelevant due to causality, 
as stressed in section IB (see also hgure 3). However one expects the Langevin equation (Ill) to describe a non 
Markovian process in which the memory kernel persists for a duration of the order of Tq,. How are we going to ignore 
the times t < to if the kernel M extends to the remote past? 

An equation like (111) may be thought of as having originated in a system coupled linearly to a bath of harmonic 
oscillators, a la Zwanzig [31, 63]. Let us consider the Hamiltonian evolution of this system, described by coordinates 
denoted collectively by F = {/i,pa,9a}, according to the Hamiltonian^®: 


-f?tot(r) — 2 ^ + V{^ + X) + 


liq 






2mr, 


Qa 


Helliu) 


HsiT) 


( 112 ) 


Note that sdetQ = numerical constant at the saddle point level 
if the system is causal. 


If we restore the inertia term in the dynamics, we would add to 
r one more coordinate momentum of the particle of position 
Similarly, Hy would have one more term p^/2m. 
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{cq, Wq, ma} are chosen suitably to reproduce dissipation and noise terms in (111). As before, assume we start in the 
remote past tp at a point in phase space Fp and let us distinguish two times to and ti, such that ti > to ^ tp. We 
can rewrite averages over the effective Langevin process as averages over the Markovian process T{t), by defining the 
transition probabilities, for fixed Fp and Fi, 




f’y(ri,ti|Fp,tp) = 


D[F] pv [F(t)|Fp, tpj 


(113) 


where pv is the transition probability density of the whole system {effective particle} U {bath}. The marginal of 
/9v[F(t)|Fp, tp] over all trajectories of the bath degrees of freedom {pa,<la} is the marginal of py[/x,/i|/Xp, tp] over 
trajectories of the response field p of the MSRDDJ probability density in (108), namely in compact SUSY notation 


with fixed initial conditions at tpi 

I pti pti 


(114) 


(a.6Mb)-h}dayMa)+A) 

We can now write, using Markovianity: 

f’y(ri,ti|Fp,tp) = J dFo Py(Fi, ti|Fo, to)Py(Fo, to|Fp, tp) 
where^° dF = IlaFor to — tp ^ Ta the system relaxes, i.e. we may assume that Py(Fg, to|Fp, tp) = 


(115) 


Fy^(Fo), the equilibrium distribution in the phase space of the whole system: 


Py'^iTo) = 


c-/3iXtot(ro) 


= p^^{po)PS^iro) 


ZvZb 

with the marginals describing respectively the effective particle 

g—Mo/2'^iiq —/3 F(aio+A) 


Py'^M = 


Zv 


Zy 


and the bath 


^B"(ro) = 


„-/3ffB(ro) 


with 


,-PHuiro) _ 


n 


27T 

PuJa 


We may now get rid of the bath degrees of freedom as in Zwanzig’s calculation, here integrating them out: 

Py(Fi,ti|Fp,tp)= [ ^PPipo) fU^pMP^^iTo) B[T]py[Tit)\To,to] 

J v27r J ^ Jro,to 


(116) 


(117) 


(118) 


(119) 


Following Hanggi [31], the last term (underbraced) can be interpreted as the transition probability of a Langevin 
process, with full equilibrium at to between the harmonic oscillators and the effective particle, given by Hb, i-e. not 
neglecting the coupling to the particle. This generating functional reads, using again MSRDDJ: 


l = y dTiPy{Ti,ti\Tp,tp) = J J B[p,p,]py[p,p,\po,to] 

where now the effective Langevin process ruling the dynamics is: 

7/i(t) = -^p(t) - P [ dt' M{t- t')p{t') + F(p(t) + A) + C(t) 

^liq Jtn 


( 120 ) 


with ^^(t)} = 0 and (^C,{t)(^{t')') = 2'^T5{t — t')M{t — t') 
p{to) = po is picked with the equilibrium measure P^{po) 


( 121 ) 


Notice the lower limit of the friction kernel and the fact that there is no extra term, called ‘initial slip’ in [31], due to 
the somewhat unusual ‘conditional average’ over the perturbed bath, given by Hb- 


Constants are arbitrary here, they can be absorbed in the nor- 
malizations. 
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FIG. 3. The procedure can be more intuitively interpreted as a resummation before to of all possible trajectories starting from 
a fixed initial position at tp (thanks to the relaxation towards equilibrium) and after ti (owing to causality), leaving us with 
finite times motions in the interval [to,ti]. 


E. Relaxation at long times 


From now on we will be able to compute physically relevant observables, and with this in mind we will use standard 
MCT methodology [1, 2], At long time difference t — t' —>■ oo, equation (106) gives 


0 = -^C{oo) - l3M{oo)[C{oo) - R^] 


liq 


Coo) = , , « l3-M(oo) = ‘ 


1 + Aliq/32M(oo) 


Aliq Aliq - (7(00) 


( 122 ) 


In the long time limit, decorrelation occurs if we assume ergodicity (which will be questioned at the dynamical 
transition, see V A), the average in (110) splits and by TTI we have to compute equilibrium averages: 


M(oo) = I y dAe^-^>-/2 + A))' 


(123) 


Using the procedure of the last subsection, {F{g,o + A)) = J P^{p,o)F{fj,o + A). Note that in equilibrium we can 
show through an integration by parts that, for any observable 0{yLo), 


^(5|) = £;<'‘»°>-7ta+A)o) 


In particular by choosing 0 = 1 we obtain 


{F{po + A)) — —— (^o) 

^liq 


For simplicity we focus here on hard spheres. For this potential the normalization reads: 

= yA)“ 0 o(-A/y^) with 00 (x) = ^ ^ erf(-x/v^) 
J-\ V27r ^ 


Given that 


(mo) = 




y /Aliq 0o( — A/ ^ Aliq) 

after a short computation we find in the limit i? —>■ oo: 


A 


liq 


„-AV2Aii< 


2^ 0o(-A/v^) 




ip e 


-Aic/4 


Aliq-foo di/fr ^Aliq 


(124) 


(125) 


(126) 


(127) 


(128) 


which shows with (122) that both M{oo) and C{oo) go exponentially to zero for A^q —>■ oo, as one would expect. 



















26 


F. Getting rid of the sphere S‘‘(R): infinite radius limit 

Applying the method in IV D, the equation for M (110) reads 

- t') = 11 dA 1 ^ (Fifiit) + X)F{fi{t') + A)) (129) 

where the average is computed over the Langevin process (121). Fixing A, we make the change of variables^^ h{t) = 
fi{t) + A centered at the wall. The normalization of is unchanged and 


Py'^iho) = 


,-{ho-\f/2Aii^-l3V(ho) 


V V-G/ 

We can simplify the expression of M with a saddle point method^^ for A^q —oo, setting a = A/Auq-. 
M(t-t') = {F{h{t))F{h{t'))) 


(130) 



(131) 


since a = 1 at the saddle point. We replaced the normalization by 

= 


^1*0 „-/Xo/2Aii<,-/3y(/xo+QAH,) 




Aliq—>-00 


a/A 


liq 


(132) 


since the potential goes to zero at long distances. The Langevin equation is affected by the change of variables only 
with an additional term Ta ~ T: 


(133) 


Our problem is mapped onto a one-dimensional diffusion with colored noise (as usual in mean field [22]) and a 
harmonic effective potential w{h) perturbed by the spheres’ repulsion (cf. hgure 4). The underbraced terms -the 
harmonic potential well- are negligible for finite times, but necessary to confine the system: they represent the ‘box’. 




FIG. 4. Effective potential landscape. If V is hard, there is an infinite wall at h = 0 prohibiting any motion in the h < 0 half 
line. If the constraint is softer, as drawn here, motion is possible for h < 0 but is rather unlikely. 


Note that h is called y in the main text. 

22 From the Langevin equation, the bracketed term 
(^F(h(t))F(h(t'))') depends upon a but it cannot be expo¬ 


nential in Aiiq and as a consequence gives no contribution to 
the saddle-point equation. 
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G. The ‘Lagrange multiplier’ 


Plugging in (105) the value at equal times^^ M{Q), which can be computed at the equilibrium using TTI, we have 


- ^ +/3M(0) = I / +/ 3 I / d/io 

^liq ^ J — 00 ^ J — 00 


(134) 


since V(ho) = 0 for ho > 0. We will focus in the following on hard spheres and assume a regularization of the potential 
which is of class C^, e.g. soft spheres^^ foss(^) = Kr^&{—h). We have, with integration by parts, due to the continuity 
of V and F in 0, 


T 


An 


liq 




2 


eho-hV{ho) (^F{ho) - T) 


dho e'‘«-^^(^'>)F(/io) (1 + l3F{ho)) +P dho 


+ T d/io 


(135) 


The last integral being zero for hard spheres, we conclude that the Lagrange multiplier^^ is (up to exponentially small 
corrections in A^q —)• 00 due to M(0) and Sv): 


pv 


2 


1 

Aliq 


(136) 


V. PHYSICAL CONSEQUENCES OF THE DYNAMICAL EQUATIONS 
A. Plateau and dynamical transition 

1. Metastable glassy states: plateau value 


We now look for a plateau in the dynamics: we assume a strong separation between a fast and a slow motion. We 
can split the correlations (and similarly the memory kernel) into a vibrational short-lived contribution and a slowly 
decaying function: 

C{t-t') = Cf{t-t') + C‘^{t-t'), + (137) 


each decaying on timescales Tf Tg, respectively. Let us look at intermediate times Tf t — t' Tg- The slowly 
varying functions are approximately constant at this scale, equal to the plateau value noted Cea (respectively Mea)- 
We have Cea = Auq — Aea where Aea is the plateau of the MSD. Similarly to (122), we get from (106) the relation 
between Mea and Aea: 


/3^Mea 


1 Cea 
A liq Aliq — C'ea 


1 1 
Aea Aliq 


(138) 


From (137), we can consider the Langevin noise as the sum of two independent centered Gaussian noises, a slowly 
varying one C and a fast one . In this limit, the Langevin equation (133) reads, using (138), 


jh{t) = s — ^ h{t) — P [ dt' M^{t — t')h{t') + F(h(t)) -b (t) 

Aea Jto 

with s = C + pM^Aho + T , (C'^(t)C'^(i')) = - t') + M^(t - t '), (139) 

{C,{t)C,{t')) = M‘‘{t - t') and for T/< t - t'<C Ts, ~ Mea 


23 M(0) diverges in the hard-sphere limit, which is natural given 24 actually valid for a linear regularization Vo^(h) = 

its interpretation as a force-force correlation. tii ,c 3 i n 7 i.r/zv\ 3/ 

—KhKy{—h) too. indeed we get ou = and pM(0) = 2 ^'^ ~ 

T) + 0{1/k). 

As noted in ID, u is not the actual Lagrange multiplier, but is 
related to it. 
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Following [32, Sec. 4.], s acts as a quasistatic field: for times t — t' Tg, (C^o) or equivalently s can be considered 
as quenched variables, picked with probability F'siow(s)- For t — t' ^ Tf, the process relaxes to an ‘equilibrium’ state 
selected by s, which is the actual metastable glassy state, with probability^® 


Pi{h\s) = 

Fslow('S) = 


g-/3Hi(h,s) 


Zi{s) 

' dfep 

dh. 


with Hi{h, s) = 


T 


2Aea 


— sh + V{h) 


dC (^o) - C - /3Mea/io - T ) 


IL g-Pw{ho)^ 


•v/SttMea 

“ea (^ho-- i^(s-T)^ Aii,/2-A+/io(a-l)-Aiiq(a-l)V2 


J \/^ V^SttMea 

with a = A/Aiiq. Taking Aiiq —>■ oo as in (131) provides the plateau value: 


Zv 




As an example, we restrict ourselves to the hard-sphere case. Then, with an integration by parts, 

T 


{Fih)),= J ^Pi{h\s)F{h) = — 




with Zi(s) = 


(s) = j ^ ^fSsh-h^l2A^^-mh) = 0^ e- 


setting u = VAea(/?s — 1 ) finally gives 


1 

1" du e 

/2 -(u+VA^)V2 

taeI 

00 1 

(—M — V Aea) 

1 


(140) 


(141) 


(142) 


(143) 


2. Equivalence with the one-step replica-symmetry-breaking result 


In [12, 30], the equation for the plateau value was obtained, maximizing the hard spheres’ replicated entropy in 
d —>■ oo in the glass phase where the IRSB ansatz is stable. It reads 

—^— = 0 (A) with GmiA) = 1 — m [ ,— (v — \/2^\ (144) 

A dA J VStt V / 

where 2A = Aea- Here we are interested in the dynamical transition where the replica symmetric solution is valid 
(the liquid state is our ‘paramagnetic’ state), which can be recovered in the limit m —)■ 1 in Monasson’s scheme. Then 


GmiA) = (l-m) 

m—^1 


i-k 


du 


ln0, 




-I- 0 ((m - I)^) 


Deriving with respect to A, we get 



which is the same as (143) with u = —v. 


(145) 


(146) 


One has to be careful here with the limit Aiiq oo: we can¬ 
not use directly equation (131), this is why we have to compute 
^siow oil sphere and take the infinite radius limit. In prin¬ 


ciple we would also do it for Pi(h\s) but it is subdominant (as 
in footnote 22), so that we can compute it directly with (139), 
where A^q oo has already been taken. Note that one recovers 
equation (16) of the main text by setting Y = s — T. 




29 


3. Dynamical transition 


Plotting the function ip versus Aea using (143) we get a minimum at the critical packing fraction 


~ 4.80678^ 


(147) 


where the plateau is AEA(<?d) — 1-15336. This packing fraction is larger than the best known lower bound for the 
existence of sphere packings ^ ^ 6 /e [53]. As mentioned in the abstract, this means that a hard sphere system may 
be prepared in equilibrium up to these densities in times that do not scale with the size of the box. In other words, 
packings as good as this are easy to obtain, and we conclude that this would be a constructive improvement on the 
best bound known in high d. It would require the present derivation to be turned into a rigorous proof, which seems 
feasible along the lines of [55] since our calculation, though a bit tedious, is quite elementary. 



FIG. 5. Plateau value and critical packing fraction 


B. Diffusion at long times 


From (107) we obtain an expression for the diffusion coefficient for times larger than the relaxation time but still 
A <C Aiiq. In this regime, the mode-coupling-like equation for the MSD reduces to 


'yD{t — t') = 2dT — 2d^/3 J dv M(t — v)D{v — t') 


(148) 


where D = A/d is the non-rescaled MSD, that is, the MSD of the original system of particles. Using Laplace transform 
we have 


D(p) 


1 2dT 1 2dT 

j + 2d?l3M(p) ^ + 2d^l3M{0) 


(149) 


By definition M(0) = dtM{t). A Tauberian theorem then gives the long-time diffusive behaviour of the MSD 
from the small p behaviour of its Laplace transform [64]: 


T 

D(f) ~ 2dDt with D = --- ^ roo ~ 

7 + 2d2/3/“M 


(150) 


At low density M ~ 0 and we recover the usual diffusion coefficient D = Toi the free dynamics. Upon increasing 
density, M increases and the diffusion coefficient decreases. At the dynamical transition, M displays a persistence 
of a plateau, the relaxation time Tq, oc M diverges and the diffusion coefficient vanishes. One usually defines an 
exponent 7 such that, for p —>■ 



and D 



(151) 


7 is one of the so-called MCT exponents [1, 2]. 
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C. Relation to the standard density formulation of MCT 

1. Intermediate scattering functions 

In its standard formulation, MCT provides equations for density correlators between time t and the origin = 
/ (|/9gP), where (•) is a canonical average over initial conditions, and pq = Fourier 

transform of the particle density [1, 2], This correlator thus reads 






= 1 


(152) 


We can define, as in IC 1 and similarly to standard liquid theory [19], the (non-averaged) local densities of trajectories 
p^'^^[x,x] ='^S{x - x,)S{x - Xi) , p^‘^'>[x,y,x,y] ='^6{x - x,)S{x - Xi)6{y - Xj)6{y - Xj) (153) 

i i^j 


SO that the intermediate scattering functions can be written as 


+<(' 5(0 « (^J B[x,x]p 




+ 


yD[a;,y,x,y]p(^)[a;,?/,x,y]e*«l®(‘^ (154) 


where the self {i = j) and distinct {i 7 ^ j) parts are defined. Both parts can then be expressed as a function of A(t) 
through the saddle point evaluation of the integrals. In the following for simplicity we discuss only the self part. 


2. The self part in infinite dimension 


In d —)■ 00 , the self part is simply expressed in terms of the MSD. Using rotation invariance of to average over 
d-dimensional random rotations TZ: 


€i{t) ^ ( / D[a:,d]p^ ^[a;,i]e*«' 


(i)u T\„iq-[x{t)-x(o)] \ _ 


B[x, x]dn [TZx, 


cx D[x,x]p^^^[x,x] d9 


(155) 


where 9 denotes the angle between q and TZ[x{t) — a;(0)] and we used hyperspherical coordinates^^. Now we can 
proceed as in IIC and express the dynamical variables in terms of the MSD: 




(156) 


The last integral can be evaluated through a saddle-point method in d —)■ 00 . Provided |g|\/A(t)/d^/^ ^ 1, the 
g-dependent term is irrelevant for the saddle-point evaluation^® on Q. The saddle-point value of 9 is imposed at the 
equator 7r/2, 


Jo e=7r/2+e Jq 

Jm. 


d—>-oo 


(157) 


In d dimensions, hyperspherical coordinates are de¬ 
fined by = rcos^i, x'^ = rsin0icos02? ••• 5 

= r sin 01 • • • sin cos with 0^—1 ^ [0?27r[, 

The measure is — 

^d-l gijid-2 sin*^"^ 02 ■ ■ • sm9(i_2 dr d9f_i. 

If g = 0(d3/2), ^hich, a priori, does not represent any physical 
distance (much less than the typical zone spanned by vibrations 


inside a cage q = 0{d)), the saddle-point values of 0 are dif¬ 
ferent, although also simple; however the ^-dependent term now 
contributes to the saddle-point equation on A, changing its value 
and the formula (158) does not hold anymore. 
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The remaining integral over Q and u is dealt with as in IIC and is normalized to 1, since the saddle-point is not 
affected by the last term in (156) as long as q^A{t)/d? ^ 1. Together with the normalization (j)g{0) = 1, we finally 
conclude for all wavevectors satisfying the latter condition, 


(plit) = exp 

^ d—>oo 



(158) 


First, we note that the self correlator is Gaussian in c? —)■ oo, in contrast to what is found in [15, 16]. In these articles, 
the MCT equations for the plateau value (the so-called Debye-Waller factor or non-ergodic parameter, which reads 
here (jf^ ^EA/ 2 d ^ golved numerically for the hard spheres system up to d = 800, and its shape is found 

to be non-Gaussian. 

This expression is also exact for any dimension both in the free-particle regime (lengths and time small compared to 
mean free path and collision time respectively) and hydrodynamic limit (lengths and time large compared to mean 
free path and collision time respectively) [19]. 

(158) implies, by substitntion in (107), equations for the <pg, with noticeable qualitative differences with respect to 
MGT equations (such as a non-local memory kernel M). 


3. The factorization property 

A crucial outcome of MGT is the so-called factorization property [1, 2], which allows to get MCT scaling laws. It 
states that, in the /?-relaxation window [i.e. close to the plateau), the difference between the value of the intermediate 
scattering functions and their value at the plateau can be factorized into a product of a function of the wavector only 
and a function of time only: 


cs H{q)G{t) (159) 

This property is a stringent test of MCT in simulations [27]. In d —>■ oo, the self intermediate scattering functions for 
all wavevectors are governed by a single quantity, the MSD. Close to the plateau, 5/S.{t) = A{t) — Aea is small and 
from equation (158), 




2 d 2 


SA{t) 


(160) 


We conclude that the factorization property holds in the infinite d limit. Besides, equation (158) is more general since 
it provides all orders in dA and is valid even far from the plateau. 


D. MCT exponents 


Starting from equations (106),(107), one can compute the different MCT exponents related to the approach to 
the plateau Aea — A{t) ~ or the departure from it A{t) — Aea ~ by expanding around the plateau value 
Aea [1“3]. We do not report here the full dynamical computation [39-42]. One finds that the MCT exponents are 
controlled by the exponent parameter A through the relations [1-3, 39]: 


r(i-a)2 r(i + b)2 =1 + 1 

r(l-2a) r(l-k26) ’ ^ 2 a 26 


(161) 


For hard spheres, we obtain A ~ 0.70698 which implies a ~ 0.324016, 6 ~ 0.629148 and 7 ~ 2.33786 [43]. We 
emphasize that the value of 7 is consistent with numerical results obtained in [44]. 


E. Connections with the microscopic model 

1. Correlation, response and mean-square displacement 

We wish to establish a connection between microscopic quantities and their counterpart at the saddle-point level 
in d —)■ 00 . Let us look at the MSD, but a similar reasoning can be done for correlations and responses (which are 
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related to it). Let us look at the microscopic MSD for a particle i 


1 / ^ 

= — {^{xi{t) - x^{t')y 


5{\YiZN[h]/N[h]) 


Zn 


Sh{t, t') 


(162) 


h—0 


since Zjsi[h = 0] = 1, with Zf 4 [h] (or .=.[h]) obtained by adding to the dynamical action a coupling to a generating 
field h{t, t'), 


N 

A[{xi, Xi}, h] = A[{xi, Xi}] “ X! / - Xi{t')y h{t, t') 

i—1 

This is a single-particle term, which amounts to shift the kinetic term in the following way: 


(163) 


(164) 


Then we use the above-derived d —>■ oo limit of InZ^/N ~ ln.=,/A^ = S with 5[h] = 5iG[h] -I- 5int, where only 

N—^OO 

the ideal gas term depends on h: 5ig = — / D[a;, a;]p[a;, i:]($[a;, x, h] -|-lnp[x, x]). Consequently, going to generalized 
spherical coordinates. 


(Xi(t) - Xi(t'))' 


SSioih] 


Zn Sh{t,t') 


(i$(Q"P,/i) 


h=o 5h(tA'') 


h—Q 


6h{t, t') 


/ dudu' D{u,u')h{u,u') = 

J h=0 


We deduce that the adimensional rescaled correlation functions 

^x,(t)-x,(t') ~ 


a-^N ^ ' d^oo a'^N ^ 5W(t') d^oo a'^N 

i i,fj. 2 V / I 

(with external fields hi) are non-fluctuating, imposed by their saddle-point value. 


(165) 


V|x,(t)-x,(t')P ~ (166) 


2. Force-force correlation and its relation to the memory kernel 

Here we establish a connection between a microscopic force-force correlation (j2i<j , where = 

—Vy(xi — Xj), and the memory kernel M. To generate such terms we will use a random shift similar to the MK 
model [29], and once again resort to the SUSY notation for compactness. We will thus consider a shift of vector 
Aijg{a) for each pair of particles, where Aij are Gaussian centered random vectors in d dimensions, of variance T,\, 
independent and identically distributed. We will note VA = /\/2^j^ their common measure, 

g is a scalar time-dependent external field that will be sent to zero in the end, in order to recover the original model. 
Again for compactness, we will note averages over the Aij by an overbar, i.e. AC = 0 and A^jAf.^ = 

The dynamical action becomes 


N 1,N 

A[{a:J, {Ag},^] = ^ d>(ic*) 4-^ / daV{xi{a) - Xj{a) d- Aijg{a)) (167) 

i=i Kj 

We still note Fy{a) = —V^U(a;i(a) — Xj{a) d- Aijg{a)), knowing that we recover the previously defined force by 
= FC(a) 1^ where the 0 stands for g and all Grassmann variables being sent to zero. Let us compute the second 
derivative of the generating dynamic functional ZN[g]- 


hZ^jg] 

6g{a) 


i<j 


5^Zi,[g] 

5g{a)g{b) 





^ A(;.A^,F^(a)F,-,(&) + d(a, b) AAA\XFtya) 

i<j i<j 

k<l 


(168) 
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where repeated Greek indices are summed over. Sending the external field to zero and averaging over the random 
shifts directly give 5ZM[g]/5g{a) = 0 and 


6"^ZN[g] 
5g{a)g{b) ^ 



Zn 


(169) 


which is the original force-force correlation looked for. As in [29], one can compute the average Zi^\g] introducing an 
averaged Mayer function 


/y[s] 


= j'DA 


Xj{a) + Ag{a)) 



g-/ daV(xi{a)-Xj{a)+Ag(a)) 


1 


(170) 


For a non-zero g{a) and large enough T,a, we still have the same crucial fact that fij[g] = 1-1- 0(Vd(r)/V) where F 
is a typical length of a trajectory (for finite times), due to the requirement that two trajectories overlap to feel the 
effect of the potential. As a consequence, we can repeat the MK computation and obtain in the thermodynamic limit 
ZN[g] = with the action 

S = - J Dx p{x){\np{x) +‘^{x)) + ^ J D[x,y]p{x)p{y) J VA f{x - y + Ag{a)) (171) 

The difference here is that we cannot simplify further by translation invariance since the shift is time-dependent. 
Nevertheless, we can still compute derivatives of reminding that due to probability conservation Zj^lg] = 1: 


5Zn 

5g{a) 


0 


= N 


dS 

6g{a) 


0 


= 0 


5^Zm 

-N 

+ N^ 

SS 

-N 

5g{a)g{b) 

0 bg{a)g{b) 

0 Sgia) 

o^9{b) 

0 6g{a)g{b) 


(172) 


We will use, as in IID, the d —> oo limit to compute S. We know in this limit that the trajectory density p{x) 
is determined, at leading order, by a saddle-point equation implying only the Jacobian of the change to generalized 
spherical coordinates, hence it is independent of the external field at this order. Thus we can overlook this dependence 
and compute derivatives of the Mayer function, leading to, after averaging: 


6^S 

6g{a)g{b) 


= 


N 


g=0 


J D[x,y]p{x)p{y)e y)A) ^p(^x - y){a) ■ F{x - y){b) -f 5{a,h)V ■ F{x - y)(a)] 

(173) 


We are interested in the “boson-boson” part, so we can focus on the first term only, which reads 


F{x - y){a) ■ F{x - y){b)L = 


\x - y\{t)\x - y\{t') 


— V (|a; - y\{t)) V' (|a; - yl(t')) 


(174) 


The content of subsection IIB is that the trajectory of (x — y)(t) ~ X plus a correction that is in l/d. For any finite 
time, the vector X can be considered to be constant, with |A| = cr and on average (A^)^ ~ a^/d. This is because 
particles do not move by 0(1) in a finite time. The correction term has zero average and \(x — y){t)\ = ct(1 J- h(t)/d). 
Recalling the definitions v\h) = V{a{l + h/d)) and F{h) = —V'{h), we have F{h) = + h/d))a/d. Using 

this and the same analysis^® of this two-body term as in IID, we get, at leading order in d —)■ oo: 


S'^ZNig] 


5g{a)g{b) 


= 


Nd^p 


J dhoe-^^^'^A {h{t)) F {hit'))) = El^M{t,t') 


We conclude from (169) 




Zn 


i<j 


(175) 


(176) 


The normalization, giving a factor proportional to the packing 
fraction as in IID 2, follows from a similar analysis since, though 
the Mayer function f(x — y) constraining the trajectories to be 


close is not there anymore (it is replaced by 1 + /), F(x — y) has 
the same role. 
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We make a final comment about two-body correlations that have the same structure as the interaction term. The 
truncated virial expansion in c? —)■ oo tells us that the two-body density of trajectories, which is the average of 
defined in (153), is simply given by^*^ [65, 66] 

p^'^'>[x,y,x,y\ = (p^‘^\x,y,x,y\\ = p[x, x\p[y, y] {l + f[x - y,x - y\) = fV 2 p(a 7 )p(y)e" (^ 77 ) 

\ / Zn 

We obtain the relation, for a function O, 

^ X! (a;*j(0))^ ^ ^ (x(f) - yit)) O (a:(0) - y(0)) 

~ iV y D[x, x]D[y, y]p[x, x]p[y, y] (l + f[x-y,x- y]) O {x{t) - y{t)) O (a;(0) - y(0)) 

(178) 


To this kind of function we can apply the same reasoning as in subsection IID, with the replacement 
/ —> (1 -I- f)O{t)O{0), if O rejects trajectories that do not get close enough^^, as was the Mayer function’s role, 
now replaced by 1 -I- / = e~^ (which is 1 for most trajectories), giving an average over the effective dynamics with 
potential. This argument gives more directly Eq. (176). However one has to be careful for correlations of a higher 
number of particles, as in the next subsection. 


3. Stress-stress eorrelation and its relation to the memory kernel 


Here we repeat the former procedure to obtain the link between the correlation of the off-diagonal (p ^ v) compo¬ 
nents of the stress tensor, which reads [19, 45] 


(J ^1/ — ^ ^ i^Xj Xj)^V Vi^Xi 

i<j 


(179) 


Once again, we use external random fields to generate the correlation function. This amounts to turn the dynamical 
action into 


N 

A[{x,},{A,j,B,j},gi,g2,hi,h2] =^$(x*) - 

Z =1 


i<j ' 




f r 

2 E (aji - Xj Aijhi) (a) -^V {x^- Xj -f (a) 


(180) 


where gi,g 2 ,hi,h 2 are external d-dimensional fields and {Aij,Bij} are one-dimensional independent identically dis¬ 
tributed centered Gaussian random variables of variance and respectively. Note that when the external fields 
are zero, we recover the original model. Using the shorthand notation 9i234 = 5^/5gi{a)g2{b)h'{{a)h2{b), we have, 
microscopically, 

_ r ^ I _ 

9i234^Ar|gi=h.=o = / Daji e ^ -Aij{x^-XjY{a)Bki{xk-xiY{b)A^nFliin{a)BpqF^q(b) 

i=l i<j 


i<3 

k<l 

m<.n 

p<q 


y2 y2 


/^(aji - XjY{a)ixk - xiY{b)Fl^Y(^)F^i{b)\ = (cr^^(a)cr^^(6))^^ 

\ i<j / 

k<l 


( 181 ) 


The N'^ factor only comes from the choice of a different normal¬ 
ization of p[x, x] its definition (29), with respect to the one used 


in liquid theory. 

This is important for e.g. the normalization as in IID 2. 
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InZ^v generates the connected correlation functions (cumulants), so that, using another shorthand notation 9" relative 
to any derivative with respect to the fields involved in 9i234: 


9i 234 In = 9^ In = == Zjq — n ZmO^Zm — ^ 

1 term ^ 4 terms ^ 


O^ZnO^Zm 


3 terms 


Zn 


O^ZnO^ZnO^Zn 


6 terms 


6 


(182) 


Zn 


ZnQ^ ZnO^ ZnO^ Zn 


1 term 


By isotropy (or average over the disorder), when evaluated at zero external field, all terms containing a first derivative 
are zero. Furthermore, second derivative terms are also zero, either due to AB — 0 for terms involving different times, 
or by isotropy for terms involving different indices^^ ^ ^ v. We are thus left with only one term, which is 9i234.^Ar 
since Zn = 1- Once again as in [29], one can compute In ZAr[gi, ^ 2 , hi, h 2 ] introducing an averaged Mayer function 


= / VAVB f{xi{a) - Xj{a),gi,g2,hi,h2) 


= / VAVB 


da (4lgi (a)+i3g2(a))'(xi —Xj )(o)—i V(xi —a;j+Ahi)(a)—i V(xi—Xj+Bh. 2 )(o) ^ 


(183) 


Note that, for finite times, e~^ /da v(xi-xj+Ah,i)(a) _ ^ js 2 ero except if hi is for some time approximately in the same 
direction as Xi — Xj, and in that case it is of order 0(T/L) where once again F is the typical length of a trajectory 
and L'^ ~ V, for large enough This implies, for fixed non-zero hi and small gi = egilZiA with gi of order one, 


y ~ y 2?^g/do4lgi(a).(xi-Xj)(o) 

= gS)i[/da4lgi(o).(xi-X3)(a)]^/2 (184) 

= l + 0{e^)+0{T/L) 

One must be careful with the order of limits, the reasoning is the following here: 

1 . we fix hi (non-zero) 

2 . we fix gi = egijYiA with gi of order one 

3. we take large so that Ahi covers the whole line spanned by hi, and T>A ~ dA/L 

4. we then take e —>■ 0. In the end we will set all external fields to zero so defining them only in the neighborhood 
of 0 is enough^^. 

The same holds for the terms involving g 2 and h 2 , with B instead of A. This way, = 0(£^) -I- 0(r/L) is small and 
we may repeat the same steps as in [29], so that InZ^v = NS, where 

S = - j B)x p{x){\Ap{x) + ^{x)) + ^ j B)[x,y\p{x)p{y) J VAVB f{x - y,gi,g 2 ,hi,h 2 ) (185) 

We only need to compute, to leading order^"^ 

di234S\g.^h.=o = y J D[x,y] p{x)p{y)e~ - yY{a)B{x - yY{h)AF''{x - y){a)BF''{x - y){b) 

= y D[a:,y] p{x)p{y)e~ ^ - yY{a){x - yY{b)F''{x - y){a)F''{x - y){b) 

(186) 


These are averages of an expression proportional to {x — y)^(x — 
y)^ which is its own opposite when rotating by an angle 7r/2 in 
the (^, u) plane. 


We can choose e = 1 /Sa if do not wish to introduce an additional 
parameter. 

As in VE2, for d ^ oo we can ignore the external fields depen¬ 
dence of the trajectory density p. 
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Likewise (174), we have for d —)■ oo 


f — v)^ {t'] (x — 

{x - yY{a){x - yY{b)F''{x - y){a)F^{x - y){b)\g = ^“ 2/1W) (1^ “ 2/K^')) 

ixYYxn\. 


|X|2 

F (hit)) F {h{t')) 


V' {\x-y\it))V' {\x-y\Y)) 


(187) 


We deduce, as in the last subsection. 


5^ \nZN 


= jdhoe-^'"^^«'^(^F{h{t))F{h{t'))'^=^=^^dNM{t,t') (188) 


^9i{a)92{b)K{o)K{^) Ig 

and from (181) and (182) 


M{t,t ) — 


( 189 ) 


We conclude that the memory function coincides with the force-force and stress-stress correlations. 
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